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1.  Introduction 


This  is  the  second  in  a  sequence  of  reports  that  present  results  of  studies  of 
electromagnetic  properties  of  finite,  three-dimensional  spherical  systems  1 1 |*.  Al¬ 
though  this  particular  geometry  is  an  extremely  simple  idealization,  many  of  the 
results  obtained  in  the  course  of  the  research  are  quite  general.  Indeed,  they  pro¬ 
vide  considerable  insight  into  the  electromagnetic  characteristics  of  real  systems 
of  much  more  complicated  geometry  and  structure  which  do  not  readily  lend 
themselves  to  analytic  resolution. 

In  the  study  presented  here,  we  investigated  one  important  effect  of  concentrically 
enclosing  a  spherical  lossless  dielectric  within  an  ideal  conducting  spherical 
cavity.  This  system  is  clearly  of  intrinsic  importance  on  its  own.  In  addition,  it 
provides  an  opportunity  to  understand,  by  comparison  with  a  previously  reported 
study  1 1 1,  how  the  effect  of  loading  a  spherical  cavity  with  a  dielectric  sphere  dif¬ 
fers  from  imbedding  a  concentric  spherical  conductor  in  that  cavity.  Furthermore, 
the  studies  of  such  systems  of  very  simple  geometry  will  subsequently  serve  as  re¬ 
liable  reference  points  for  studies  of  systems  of  successively  increasing 
complexity. 

The  electromagnetic  property  we  focus  on  in  this  report  is  the  set  of  resonant  fre¬ 
quencies  of  the  dielectrically  loaded  conducting  spherical  cavity.  The  dependence 
of  these  physical  variables  on  size  and  dielectric  strength  will  be  presented.  In 
section  2  we  briefly  present  the  formal  development  for  the  eigenvalue  problem  of 
the  system.  In  a  previous  report  |I|  we  gave  a  more  complete  formal  develop¬ 
ment.  The  details  of  the  geometry  are  given  in  section  2  also.  Transverse-electric 
(TE)  and  transverse-magnetic  (TM)  modes  (i.e.,  relative  to  the  radial  direction) 
are  separable  and  uncoupled  for  this  system.  Section  3  contains  the  resonances 
for  the  TE  modes  of  the  system  as  a  function  of  radius  of  the  interior  dielectric. 
Effects  of  varying  the  dielectric  strength  for  the  eigenvalues  of  the  TE  modes  are 
included  in  this  section.  Corresponding  results  for  the  TM  modes  are  given  in 
section  4.  Finally,  in  section  5  we  present  a  summary  discussion  of  all  the  results. 


References  are  listed  at  end  of  report 


2.  System  Geometry  and  Formal  Developments 


System  Geometry 

Figure  1  illustrates  the  geometry  of  the  system  undergoing  investigation. 
It  consists  of  a  simple  spherical  conductor  of  radius  a  containing  a  con¬ 
centric  lossless  dielectric  sphere  of  radius  b  and  dielectric  constant  e2. 

The  remaining  volume  contains  a  lossless  dielectric  of  dielectric  constant 
£,.  The  relative  dielectric  strength  is  denoted  by  er,  where 


er  =  £2/ev 


We  shall  in  this  report  consider  only  the  case  where  the  magnetic  per¬ 
meability  is  the  same  for  both  dielectrics;  i.e.,  |i,  =  |i2. 

Quantities  related  to  the  region  between  the  outer  boundary  of  the  inner 
dielectric  sphere  and  the  conductor  will  be  denoted  by  the  index  1,  either 
as  a  subscript  or  superscript.  The  variables  related  to  the  volume  occupied 
by  the  inner  sphere  will  similarly  be  denoted  by  the  index  2. 


Figure  1.  Schematic 
of  the  cunducting 
spherical  cavity  con¬ 
taining  a  concentric 
spherical  dielectric. 
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dielectric 
load  sphere 
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'Cavity  dielectric 


Formal  Derivation  of  TE  Mode  Eigenvalue  Relation 

The  Debye  potential  for  the  TE  mode  is  of  the  form 


^(r)  = 


Ynm (0,0) %  (r)  ,  for  b<,r<a 
Y nm  (e,0)0<2)(r)  ,  for  OZrcb 


and  satisfies 


where 


(v2  ~Y?)<t>nm(r)  =0  ,  for  b  <r  <, a 


(v2  =0  >  for  O^r  <b 


Y.  3  A  =7'V^  o) 


Since  the  angular  part  is  well-known,  we  shall  only  explicitly  consider  the 
radially  dependent  part  of  equation  (2),  which  is  taken  in  the  form 

f0i,>(Yir)  =An)‘n(Yrr)  +Bi^kn(y]r)  ,  h^rZa 


0n  A3 


0„  (Y2  r)  =  *n  (Y2  r)  . 


0  <r<b 


where  we  have  introduced  the  modified  spherical  Bessel  functions  in 
equation  (5).  Note  that  the  coefficient  of  0^2)  has  been  chosen  as  unity  for 
convenience.  Applications  of  the  boundary  conditions  of  the  fields  at  r  = 
a  and  r  =  b  will  then  determine  the  coefficients  A*1*  and  and  establish 

n  n 

the  eigenvalue  relation  as  a  function  of  a,  b,  e,,  and  £2-  The  general  ex¬ 
pressions  for  the  field  components  in  terms  of  the  Debye  potential  for  the 
TE  modes  are  given  by 


Er.nm(D  =0 

E<>nm  (r)~~  rsm0  <50  [r<1Wr'] 
fe;mn(r)  =  -J-  |rd)nm(r)J  , 


(6c) 


and 


Hr*"  (r)  =  7^r^  -f  )[r*nm(r)]  . 
We;nm (f)  =  Jujl  7  3750  (r<I>«'"(r)]  • 


j  |  ^2 

^6;«m(r)  =  J^jJ  7  J5^  [r&nm  (*")]  • 
Each  of  these  field  components  is  of  the  form 

fl  (1^(0  .  b  <  r  <,a 

l  l();wn  (r)  =  1  i() \nm  (6»0)  l  n. 

,  ()<r<b 

We  then  can  separate  variables  and  rewrite  equations  (6)  as 

Er  -nm  00  =0  , 

r 4!>>  ,  b<r^a 

(r)  =  E**.(e.«»  {' 

,  0  <r<b 
ffiSJW  .  b  <  r  <a 

Efrnm 

l£&(r)  ’  0*r<b 

[  H(rn  (r)  -  b  <,r  <,a 

Hrvm{r)=Hrvm&M  ' 

{ H™{r)  ,  0  <r<b 


f Hq),  (r)  ,  b  <  r  <  a 

(r)~  Hti-.nmiQ’Q)) 

l«S(r)  ,0Zr<b 
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^6;nm(r)  =  W6;wn  (8,<J))  - 


r «2i(r)  , b<rZa 


<i (r)  ,  0  <r  <  b 


(80 


Clearly  we  need  only  consider  the  radially  dependent  factor  in  all  further 
derivations  in  this  report.  Thus  we  shall  write: 


Ee:n(r)  = 


Er-.n  (r)  =0  . 

(9a) 

=  <J>1'V)  .  b  <  r  <a 

l4;>>  =^(r)  ,  0<  r  <  b 

(9b) 

Efr.nm  (6><t>)  | 

4';i(r)  ,  b<r<,a 

,0<r<a 

(9c) 

'H{r])n(r)  ,  b  <  r  <a 

H{^n(r)  ,  0  <r<b 

(9d) 

f«Sj,(r)  ,  b<rza 

1  H$n(r)  ,Q<rcb 

(9e) 

and 


Or> 


«£(r)  ,  0<r<a 


These  can  be  further  expanded  to 

£<-;«('■)  =0  , 


£e;«('')  = 


4!n(r)  ^L’V)  -  b  <r  <a 
=  <2)  (r)  ,  0  <  r  <  <i 


(90 


(l()a> 


(10b) 
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£*-.nW  =£0;„(r)  , 


( 10c) 


777  777777) 

W&  7^7  [Yi'-‘t>»"(Y.r)]' 


W9;n(r)  =<( 


wS(r)  =  A/S  ^7  57^77  [^'-^’(x.O] 


M2  72  r  d(y*r) 


|  777  b''C(Y2'-)]' 
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and 


(10e) 


H0,n(r)  =H».n(r)  ,  (100 

where  we  have  introduced  the  notation 

^  I  v<t>(.v)|  h  |.v<t>U)|  '  . 

Now  let  us  explicitly  consider  applying  the  boundary  conditions.  At  the 
dielectric  interface,  both  the  tangential  electric  field  and  the  tangential 
magnetic  field  are  continuous  across  at  r  =  b.  From  equations  (9b)  and 
(e),  we  obtain  then  the  pair  of  relations 


(11) 


K  \y\hin(y\t>)\'  +  B'n,\yxbkn(y]b)\'=\y2bin(y2h)]'  .  (12) 


This  pair  of  equations  is  readily  solved  to  give 

(!)__  h  ln (lib)  .  MYi^) 

Y‘  |Y  '>^(Y2^)| '  •  IYi&MYi'’)! ' 


and 


in(y\b)  .  i„(  y2b) 

\y\bin(y\h)\ '  ,  \y2bin(y2h)\  ' 


(13) 


(14) 


where  we  have  used  the  Wronskian  relation  to  simplify  the  results.  The 
coefficients  have  thus  been  determined  in  terms  of  known  quantities. 

Next  we  apply  the  boundary  condition  at  the  conducting  surface:  namely, 
the  tangential  electric  field  vanishes  at  r  =  a.  This  gives  us 

A\]]in{ Yt“)  +fl|,l)MYi<0  =()  • 

Upon  substituting  into  this  relation  from  equations  (13)  and  (14)  for  the 
coefficients  and  rearranging,  we  obtain  one  form  of  the  eigenvalue  rela* 
tion  for  the  TE  modes: 

in  (Y2 b)  rj.1  \a,b)  ~[y2b in  (y2  b)  j '  =0  .  (15) 

This  is  the  relation  whose  roots,  y^,  give  the  system  TE  eigenvalues.  The 
auxiliary  functions  introduced  in  equation  ( 15)  are  defined  as  follows: 

(a ,b)  m  in  (y, a)  kn  (ytb)~  in  (y, b)  kn  (y,a)  .  (16) 

sin(yla)\y,hkn(yib)\ '  -  kn(y,(i)\y,bin(y,b)\ '  .  (17) 

In  section  3  we  shall  continue  with  the  considerations  of  the  TE  modes. 
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2.3  Formal  Derivation  of  TM  Mode  Eigenvalue  Relation 


For  the  TM  modes,  the  Debye  potential  is  of  the  form 


(r)  = 


Y,tm (0.0) V)  .  for  b  <r  <a  , 
)'mn(0.0)V<„2)('')  ,  for  0  <  r  <  b  , 


which  satisfies  the  differential  equation: 


(v2  ~Yi)  ^nm(r)  =0  ,  for  b  <r  <a 


(V2  -  Y2)  'Fjimfr)  =  0  ,  for  0  £  r  <  a  . 


(19b) 


Once  again  we  observe  that  the  angular  part  is  well-known.  Hence  we 
only  explicitly  display  the  radially  dependent  part.  This  is 

(Vn ’(Yi r)  =Cn)/«(Yi'-)  +Dn]kn(  Yi  0  -  forb<r<,a  . 
Vn(0={  (20) 


rn  (Y2  0  =  in  (fir) 


,  for  0  <r  <b  . 


Again  the  coefficient  of  the  inner  Debye  potential  has  been  taken  as  unity 
for  convenience.  Application  of  the  boundary  conditions  for  the  fields  at 
the  conductor  and  the  dielectric  interface  will  as  before  explicitly  evaluate 
the  coefficient  Cjl1)  and  D(^  and  generate  the  eigenvalue  relation  for  the 
TM  modes. 

The  general  formal  expressions  for  the  TM  mode  field  components  in 
terms  of  the  Debye  potential  \j (r)  are 


Hr:nm(r)  =0  - 

H%\nm  ( r)  =  7^,0  ^  |''%im(r)| 


(21b) 


^$;nm(r)  ~  r  50  (rV^wn(r)l  * 


(21c) 


Ennm(r)  =j^\^2  -Y2jl^nm(r)|  , 

(2  Id) 

Ee,nm(r)  -  7  5755  |'‘H'wn(r)]  . 

(21e) 

Efcnm  (r)  ~ycoe  rsjn0  575^  l'''^nm(r>]  • 

(210 

Again  observe  that  each  of  the  field  components  in  equations  (21)  are  of 
the  form 


Ho  \nm  (r)  =1  l<)  \nm  <e,w 


fl  I')-,*'') 

ll  6,w 


b  <r  <a 
0  <r<b 


and  just  as  we  did  for  the  TE  components  we  rewrite  equations  (21)  as 
follows: 


Hnnmi  r)  =0 


He  \nm  (r)  =  He  \nm  (04»{ 


O)  ,  b  <r  <a 


{H{el(r)  ,  0 <r<b 


H$l(r)  ,  b  <r  <,a 


//<t.wn(r)=//0:nm(0,<t))^ 

{H™{r),0<r<b 


\^n(r)  ,  b  <r  <a 

Er;nm(r)  =  Er,nm  (0$)  1  n\ 

l?%(r)  ,  0  <r<b 


Eq X)n(r)  ,  b  <r  <a 


Ee  \nm  (r)  =  Ee  \nm  (0,0)' 


E%W)  -  0<  r  <  b 


(22a) 


(22b) 


(22c) 


(2  2d) 


(22e) 


£,;«('')  =  < 


E*n(r) 


and 


(r)  ~  (9)0)  ^ 


4!n(r>  ’  b<r<a 
l*£(r)  ,  0  <r<b 


Extracting  the  radially  dependent  part  we  have 

Hnn(r)  =0  , 

f  H(ll 

H&,n  (r)  =< 


<r  <a  , 


'  He-n  (r)  =V«(YiO  )  b<r<a  , 
l«£,W  =V“(»Y2'-)  ,  0 <r<b  , 

4>>  =' Vl"| -']  ,  0Sr<ft  . 

^7  [Yi  ^  Vi0  (Yi  0]  '  ,  b<r<a  , 

4£ <r>  -  ^7  [l6'V,2’<Y.  r)]  '  ,  0  <  r  <  ft  , 


/" 

■< 


^(O  =£e;„(r)  . 


(220 

HI 

(23a) 

S{r 

(23b) 

(23c) 

1 

(23d) 


(23e) 


(230 


Next  we  apply  the  boundary  conditions.  The  tangential  components  of  the 
electric  and  magnetic  fields  are  continuous  across  the  dielectric  interface 
at  r  =  b.  Then  from  equations  (23b)  and  (23d)  we  obtain,  respectively,  the 
pair  of  simultaneous  equations 


C{n)in(y\b)+D^kn(yxb)=in(y2b) 


(24) 


,v.V 

»  •  •  *- 

•  **  -V 

m 


■$ 


>'V 


wj 

vv1 
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y, /,  1^1  ^  '«(Yl  ^)]  +  Dn  y]  fj  |Yl  ^  MYl  1  <25) 


1 

ViT  Y2/> 


|Y2*<«(Y2*)r  • 


Solving  this  pair  of  equations  we  obtain,  after  some  rearranging  and  using 
the  Wronskian.  the  coefficients 


1 


y-*(  I)  _  _  * 

"  "  yih 


M Yl  h) 


El 


yibin{y2h)  .  IYi^MYi^)! 


(26) 


6.(Yl  *>)  ,  *«(Y2*>) 

Yi^‘«(Yi*)  •  ^|Y2^<«(Y2*)r 


(27) 


At  the  conducting  surface  r  =  <i,  the  tangential  components  of  the  electric- 
field  vanish.  This  gives  us  the  relation 

ClY.dMY.a)!'  +D("\y]akn(y]a)Y  =0  .  (28) 

If  we  substitute  into  this  the  expressions  for  the  coefficients  from  equa¬ 
tions  (26)  and  (27),  we  obtain,  after  some  rearranging,  one  form  of  the 
eigenvalue  relation  for  the  TM  modes  of  the  system: 

+[Yifri«(Yi*)r7i,)(M)  =0  .  (29) 


In  equation  (29)  we  introduced  the  additional  auxiliary  function 

7 4nUhb)  =|Y/<i/„(Y/<i)|  \yihkn{y,b)\  ' 
~\yibi„{yib)\  (Y/«^(Y/«)|  '  • 


The  TM  eigenvalues  were  found  by  solving  equation  (29)  for  its  roots  y,a. 
We  shall  continue  the  discussion  of  the  TM  modes  in  section  4. 


3.  The  TE  Eigenvalues 


3.1  Introduction 

Before  presenting  the  TE  eigenvalues,  we  must  first  establish  some 
conventions. 


We  introduce  the  notation 


Rsb/a  . 


Using  the  identities 


=  ~  Ryxa  , 


Yj/>  =  -  Yi<*  =/?Yi«  , 


we  rewrite  equations  (15)  through  (17)  as  follows: 
in(^fer  Ry\(>)  T nhujia)  -zrRyxa  in(yltr  ^Yi«)  A =0  , 


where 


=\]j\ajta)  =  in(yia)  kn(R  Y\  a)  -  i„  (R  y,  a)  kn  (Yi  a)  . 


Ti"  (<!,*>=  TV Wo) 

(34) 

=  *» (Yi  a)\Ry]akn (R  Yi «) | '  -  k„(y\a)[Ry\a  i„  (/?  Yt «)  ]  '  • 

Clearly,  for  fixed  values  of  er  and  R  we  can  consider  yta  as  the  variable 
in  equation  (32)  and  solve  for  those  values  of  the  variable  that  satisfy 
equation  (32).  These  values  we  call  the  "eigenvalues"  for  the  TE  modes. 
Since  the  dielectrics  within  the  cavity  are  lossless  and  the  outer  cavity 
wall  is  a  perfect  conductor,  all  the  roots  of  equation  (32)  are  purely  im- 


aginary.  Thus,  the  TE  eigenvalue  spectrum  for  each  value  of  n  will  be  a 
doubly  infinite  discrete  set  of  ordered  real  numbers,  each  of  which  we 
denote  by  k^a  where  jk^  (a)  is  a  root  of  equation  (32)  for  fixed  n  (by 
allowing  yj,  =  -fy,  and 

<k^-^  f°r  p  <p >  p  andn'=  2 . 

Further  each  value  a  is  a  function  of  the  relative  dielectric  strength, 
er,  and  of  the  ratio  of  inner  dielectric  radius  to  outer  cavity  radius,  R : 

•  <«> 

First  we  consider  two  limiting  cases: 

•  The  inner  dielectric  has  b  =  0.  In  this  limit  the  eigenvalues  are  exactly 

the  same  as  those  for  the  well-known  spherical  cavity  |2|,  filled 
homogeneously  with  dielectric  e,.  Thus,  for  each  n  the  a  are  found 
from  the  roots  of  /„  =  0. 

•  The  inner  dielectric,  e2,  fills  the  entire  cavity,  so  that  b  -  a  (or  R  =  1). 
In  this  limit  we  find  the  eigenvalues  Y/’d  by  merely  taking  those  al¬ 
ready  obtained  for  the  cavity  filled  with  dielectric  of  strength  e,  and 
dividing  by  Ve^.  That  is  to  say, 

k{"(R=0)a  =  \/y[trk{"{R=\)a. 

Clearly,  no  more  effort  is  required  to  obtain  the  TE  eigenvalues  when  the 
inner  dielectric  sphere  has  zero  radius  than  in  the  other  limit  when  the  in¬ 
ner  dielectric  has  expanded  to  fill  the  cavity  completely. 

We  note  that  if  er  >  1  (i.e.,  the  inner  dielectric  constant  exceeds  the  outer, 
surrounding,  dielectric  constant),  the  eigenvalues,  in  general,  shift 
downward  as  the  inner  sphere  fills  the  cavity. 

This  is  of  course  to  be  expected.  (However,  it  still  remains  to  be  seen 
how  the  TE  eigenvalues  behave  in  between  these  limits.)  In  the  situa¬ 
tions  where  er  <  1  (i.e.,  the  dielectric  constant  of  the  inner  sphere  is  less 
than  that  of  the  surrounding  dielectric  in  the  cavity),  the  reverse  behavior 


is  of  course  to  be  expected.  That  is,  when  the  inner  sphere  fills  the 
cavity,  the  TE  eigenvalues  y:%  take  on  the  values  that  are  well-known 
for  the  cavity  with  a  homogeneous  dielectric  in  its  interior.  As  the  inner 
sphere  shrinks  in  radius,  the  TE  eigenvalues  reduce  to  this  latter  set  mul¬ 
tiplied  by  the  factor  l/^[Tr  >  1.  Hence,  we  now  observe  that  in  the  ab¬ 
sence  of  the  inner  dielectric  sphere,  all  the  TE  resonance  frequencies  are 
shifted  to  lower  values.  (Again,  we  have  yet  to  see  the  behavior  in  the  in¬ 
termediate  region.) 

We  now  proceed  to  examine  the  behavior  of  the  TE  eigenvalues  as  a 
function  of  the  radius  of  the  concentric  dielectric  sphere  in  the  cavity  for 
several  representative  choices  of  dielectric  strength.  Though  other  repre¬ 
sentations  can  be  construed,  we  choose  to  keep  the  relative  dielectric 
strength,  e2/e ,,  fixed,  and  we  solve  equation  (32)  for  the  eigenvalues,  y,a, 
while  varying  the  normalized  inner  radius,  b/a.  The  solutions  in  this  rep¬ 
resentation  group  themselves  naturally  into  families  of  curves  in  accor¬ 
dance  with  the  second  subscript,  p.  We  shall  present  the  solutions  so 
grouped  into  these  families. 

The  TE  Eigenvalues  for  er  =  3.00 

The  first  detailed  case  we  shall  consider  is  that  for  which  the  inner 
dielectric  sphere  has  a  dielectric  constant  of  e2  =  3.00e,.  Note  that  e,  is 
assumed  to  be  the  dielectric  constant  in  the  cavity  in  the  absence  of  the 
concentric  dielectric  sphere. 

In  figure  2  the  results  obtained  from  equation  (32)  for  k^  (ft;er  =  3.00)  a 
are  shown  plotted  against  R  for  the  family  p  =  1  and  n  ranging  from  1  to 
5.  Table  1(a)  contains  numerical  values  for  these  eigenvalues.  Examina¬ 
tion  of  these  k^l(Rxr)a  curves  shows  the  following  behavior  between 
the  limiting  cases  discussed  above.  For  a  given  value  of  the  index  n,  the 
eigenvalues  at  first  decrease  in  value  slowly  as  the  normalized  inner 
sphere  radius,  R,  grows.  The  greater  the  index  n,  the  slower  this  decrease 
in  eigenvalue  becomes  and  the  larger  the  range  in  R  over  which  this  in¬ 
itial  behavior  persists.  Immediately  following  this  characteristic  initial 
behavior  with  increasing  radius  b  (or  equivalently  ft),  a  speed-up  in  the 
rate  of  change  of  k(f]Ja  occurs.  This  increase  in  the  rate  at  which  the 
eigenvalue  decreases  becomes  more  marked  as  n  increases.  There  is  a 
clearly  defined  knee  in  the  k'Jja  eigenvalue  trajectories  for  n  >  1  marking 
where  the  curve  departs  from  the  slowly  varying  to  the  quickly  varying 
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Figure  2.  TK  eigenvalues  (/?  ;e,  =3 )a  for  the  spherical  cavity  of  radius  a  filled  with  dielectric  £,  inter¬ 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  36,. 

behavior.  For  large  n  this  knee  becomes  quite  a  sharp  bend.  Following 
these  noteworthy  characteristics  the  k^a  trajectories  show  a  slowing  up 
in  the  rate  of  decrease  of  the  eigenvalues  as  the  inner  sphere  continues  to 
grow.  For  large  R  values,  corresponding  to  a  high  factor  of  filling  of  the 
cavity  by  the  inner  dielectric  sphere,  we  again  observe  a  very  slow 
change  in  the  eigenvalues  as  R  grows  further  towards  unity.  As  R 
reaches  unity,  the  eigenvalues  reach  their  limiting  values  of  1/V3  times 
their  initial  values.  Clearly,  three  distinct  regimes  are  evident  from  the 
k"l  (/?;er=3)  a  trajectories.  These  correspond  to  the  presence  of  a  small, 
intermediate,  or  large  concentric  inner  dielectric  sphere. 

The  set  of  TE  eigenvalue  trajectories,  k^j  (R  xr) ,  raises  two  questions. 
First,  why  does  the  eigenvalue  remain  nearly  constant  for  the  higher  n 
values  over  such  a  substantial  range  of  R ?  Second,  why  does  it  then 
change  so  rapidly  with  increase  in  size  of  the  dielectric  sphere  in  the  inte¬ 
rior?  Clearly  this  behavior  will  be  even  more  pronounced  for  even  higher 
values  of  n.  The  answers  to  these  questions  can  be  found  in  an  examina¬ 
tion  of  the  corresponding  spatial  distribution  of  the  Field  components 
within  the  dielectrically  loaded  cavity.  This  subject  will  be  taken  up  in 


Table  la.  Numerical  Values  of  the  TE  Eigenvalues  k^(R\Er  =  3)  a  for 
the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 

R=b/a 

n«1 

n*=2 

n=3 

n=4 

n-5 

0.05 

4.49338 

6.76346 

6.98793 

8.18256 

9.35581 

0.10 

4.49254 

5.76343 

6.98793 

8.18256 

9.35581 

0.20 

4.46441 

5.75970 

6.98749 

8.18251 

9.35581 

0.30 

4 . 27463 

5.69606 

6.97091 

8.17856 

9.35491 

0.40 

3.83541 

5.31497 

6.74779 

8.07928 

9.31677 

0.50 

3.38457 

4.66134 

5.97362 

7.30555 

8.64409 

0.60 

3.04328 

4.10949 

5.19990 

6.30710 

7.42375 

0.70 

2.81017 

3.71817 

4.63497 

5.56076 

6.49309 

0.80 

2.66886 

3.46968 

4.26308 

5.05552 

5.84923 

0.90 

2.60509 

3.34950 

4.07203 

4.78186 

5.48385 

0.95 

2.59570 

3.33051 

4.03971 

4.73248 

5.41376 

Table  lb.  Numerical  Values  of  the  TE  Eigenvalues  {Rxr  =  3)  a  for 
the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R=b/a 

n*1 

n*»2 

n*3 

n=4 

n»5 

0.05 

7.72486 

9.09501 

10.41712 

11.70491 

12.96653 

0.10 

7.71219 

9.09434 

10.41709 

11.70491 

12.96653 

0.15 

7.61736 

9.08326 

10.41601 

11.70481 

12.96652 

0.20 

7.26888 

8.99890 

10.40211 

11.70273 

12.96623 

0.25 

6.70702 

8.61623 

10.29108 

11.67880 

12.96135 

0.30 

6.25514 

7.91051 

9.74735 

11.46473 

12.90377 

0.35 

6.01878 

7.32347 

8.86154 

10.59696 

12.34333 

0.40 

5.94320 

7.01957 

8.24572 

9.64212 

11.17173 

0.45 

5.93663 

6.93209 

7.97131 

9,08310 

10.28875 

0.50 

5.90684 

6.92584 

7.92209 

8.93074 

9.96407 

0.55 

5.78584 

6.87104 

7.90644 

8.92085 

9.93401 

0.60 

5.58018 

6.69032 

7.76864 

8.82638 

9 . 87332 

0.65 

5.34200 

6.42117 

7.48287 

8.53538 

9.58448 

0.70 

5.10996 

6.13349 

7.14235 

8.14328 

9.14047 

0.80 

4.72953 

5.63503 

6.52371 

7.40152 

8.27184 

0.90 

4.50832 

5.32483 

6.11855 

6.89681 

7.66396 

0.95 

4.46713 

5.26212 

6.03067 

6.78049 

7.51630 

Table  lc.  Numerical  Values 
the  Spherical  Cavity  Loaded 


of  the  TE  Eigenvalues  ( R;er  =3 )a  for 

Internally  with  a  Dielectric  Sphere 


R*b/a 

n-1 

n*2 

n=3 

n-4 

0.05 

10.90192 

12.32290 

13.69802 

15.0397 

0.10 

10.82516 

12.31732 

13.69769 

15.03965 

0.15 

10.26354 

12.21370 

13.68510 

15.03826 

0.20 

9.43374 

11.46032 

13.48429 

15.00451 

0.25 

9.15424 

10.48578 

12.33261 

14.44791 

0.30 

9.13449 

10.21737 

11.43087 

12.94811 

0.35 

9.01077 

10.20009 

11.29319 

12.41148 

0.40 

8.61039 

9.99714 

11.24205 

12.38336 

0.45 

8.14505 

9.49576 

10.84570 

12.14966 

0.50 

7.82181 

8.99353 

10.23540 

1 1 . 50694 

0.55 

7.70064 

8.68985 

9.75206 

10.87014 

0.60 

7.69066 

8.61488 

9.54641 

10.50728 

0.65 

7 . 60805 

8 . 59006 

9.52819 

10.44599 

0.70 

7.38682 

8 . 40977 

9.39927 

10.36375 

0.75 

7.10426 

8.10593 

9.08477 

10.04735 

0.80 

6.82694 

7.78215 

8.71681 

9.63658 

0.90 

6.41015 

7.26823 

8.10480 

8.92541 

0.95 

6.31397 

7.14057 

7.94314 

8.72761 

0.99 

6.29566 

7.11489 

7 . 90888 

8.68359 

*. 


± 


R 

■=b/a 

n=1 

n=2 

n-3 

0 

.05 

14.05819 

15.51439 

16.92362 

0 

.  10 

13.74719 

15.48531 

16.92143 

0 

.  15 

12.59739 

14.86023 

16.82536 

0 

.20 

12.30726 

13.57811 

15.49753 

0 

.25 

12.23821 

13.41956 

14.58338 

0 

.30 

11.66648 

13.23988 

14.53831 

0 

.35 

11.04040 

12.48532 

14.04369 

0 

.40 

10.87308 

11.92911 

13.18152 

0 

.45 

10.82319 

11.85081 

12.84483 

0 

.50 

10.46800 

11.69153 

12.80069 

0 

.55 

9.97096 

11.19225 

12.39661 

0 

.60 

9.60650 

10.68693 

11.80285 

0 

.65 

9.49097 

10.41708 

11.38032 

0 

.70 

9.46447 

10.38642 

11.28030 

0 

.75 

9.26066 

10.24452 

11.19505 

0 

.80 

8.93346 

9.90793 

10.85978 

0 

.90 

8.32657 

9.21190 

10.07710 

0 

.95 

8.15809 

9.00530 

9.63100 

0 

.99 

8. 12147 

8.95783 

9.77147 

2‘ 


detail  in  a  forthcoming  report  giving  the  fields  in  the  cavity  for  the  TE 
eigenmodes. 

Continuing,  we  show  in  figure  3  the  results  obtained  for  the  eigenvalues 
(/?;er  =  3)ci  from  equation  (32).  Table  1(b)  contains  numerical  val¬ 
ues  for  these  eigenvalues.  Examination  of  figure  3  reveals  that  for  a 
given  index  value  n  the  initial  behavior  of  the  eigenvalues  is  similar  to 
that  of  figure  2.  That  is,  at  first,  as  R  increases,  the  eigenvalue  only  very 
slowly  decreases  from  the  initial  value  when  the  inner  sphere  is  not 
present.  This  slow  variation  persists  over  a  larger  range  of  size  of  the  in¬ 
ner  sphere  as  the  index  n  increases.  Now,  however,  for  each  value  of  n 
the  a  trajectory  begins  to  change  rapidly  with  R  noticeably  earlier 
than  the  k\^  a  trajectory.  For  the  higher  values  the  k\jl  a  trajectories  show 
a  faster  variation  with  increasing  R. 

A  new  feature  can  be  seen  in  the  curves  in  figure  3.  After  the  initial  knee 
in  the  curve,  the  eigenvalues  again  display  a  very  slow  decrease  in  value 
with  increasing  inner  sphere  radius.  As  n  increases,  this  behavior  persists 
for  a  decreasing  range  of  R ,  after  starting  at  a  larger  value  of  R  for  higher 
n.  Again,  as  R  increases  further,  we  observe  another  rapid  decrease  in 
^n2°'  which  is  more  pronounced  as  n  increases.  This  behavior  gradually 
tapers  to  a  slow  decrease  in  eigenvalue  with  R  from  about  R  =  0.8,  thus 
enabling  the  trajectories  to  go  over  to  the  known  limiting  values.  Note 
that  figure  3  clearly  indicates  two  knees  in  each  of  the  k^Ui  eigenvalue 
trajectories.  Roughly  speaking,  we  can  still  distinguish  between  the  three 
regimes — small,  intermediate,  and  large  inner  dielectric  spheres.  This 
method  of  classification  now  seems  less  clear  in  some  ways  than  earlier. 

Figure  4  shows  the  eigenvalue  trajectories  for  the  modes  k{n ]J  (Rxr  =  3)  a 
for  n  =  1,2, 3, 4.  The  corresponding  numerical  values  of  these  eigenvalues 
are  given  in  table  1(c).  Upon  examination  of  figure  4  we  observe  behav¬ 
ior  very  similar  to  what  we  previously  described  for  the  k^\a  in  figure  3. 
The  new  feature  present  in  figure  4  is  another  knee  in  each  of  the  trajec¬ 
tories  as  an  addition  to  a  behavior  similar  to  figure  3.  Once  again  we  in¬ 
dicate  that  we  shall  be  able  to  better  understand  the  structure  of  the  trajec¬ 
tories  for  the  k^\ a  eigenvalues  only  by  resorting  to  an  investigation  of 
the  distribution  of  the  field  components  for  each  eigenmode.  This  will  be 
done  in  a  subsequent  report. 
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figure  3.  TE  eigenvalues  k^  (/?  :er  =  3)  a  for  the  spherical  cavity  of  radius  a  filled  with  dielectric  £1  internally 
oaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  3£r 
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Figure  4.  TE  eigenvalue  ^  (R  ;cr  =  1)  a  fur  the  spherical  cavity  of  radius  a  filled  with  dielectric  E,  internally 
loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  3E,. 
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In  figure  5  we  present  the  final  set  of  TE  eigenvalues,  for  (R,er)  a ,  to 
be  included  in  this  report.  Table  1(d)  gives  the  corresponding  numerical 
values  for  these  eigenvalues.  Again  we  can  readily  note  from  examina¬ 
tion  of  figure  5  that  the  trajectories  for  the  k^J  a  demonstrate  the  same 
kind  of  behavior  with  R.  the  size  of  the  inner  sphere,  as  the  previous  set. 
Again  that  behavior  is  emphasized  more  strongly.  Once  again  the 
predominantly  new  feature  appearing  is  the  occurrence  of  an  additional 
knee  in  the  curves. 

At  this  stage  of  the  report,  we  can  readily  conclude  that  a  further  increase 
in  the  index  p  beyond  4  will  result  in  the  same  behavior  with  R  that  we 
have  already  seen  occurring.  Indeed  we  can,  with  reasonable  accuracy, 
predict  from  the  data  in  figures  2  through  4  the  trajectories  of  the  higher 
order  mode  eigenvalues.  With  the  aid  of  the  tables  we  can  extrapolate  to 
find  their  numerical  values  with  some  degree  of  accuracy. 


3.3  The  TE  Eigenvalues  for  er  =  10.00 

We  consider  next  the  results  of  concentrically  loading  the  conducting 
spherical  cavity  with  a  dielectric  sphere  of  large  relative  dielectric 
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Figure  5.  TF.  eigenvalue  k^l  [R  xr  =  })«  fur  the  spherical  cavity  of  radius  a  filled  with  dielectric  E,  internally 
loaded  with  a  concentric  dielectric  sphere  of  radius  ft  and  strength  E,  =  3E  . 
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constant,  namely,  e2  =  H).(X)  £,.  Figure  6  shows  results  calculated  from 
equation  (32)  for  p  =  1  and  n  ranging  from  1  to  6.  Table  2(a)  contains 
numerical  values  for  the  k^l  (R\er  =  10) a  eigenvalues  plotted  in  figure  6. 
Inspection  of  figure  6  immediately  reveals  that  the  trajectories  of  the 
eigenvalues  behave  much  like  those  shown  in  figure  2  for  p  =  1  and  e,  = 
3e(.  Each  trajectory  still  starts,  for  R  «  1,  at  the  same  value,  decreases 
slowly  at  first  with  increasing  R,  then  at  an  accelerated  rate  with  further 
increase  in  /?,  and  gradually  slows  its  rate  of  decrease  as  R  approaches 
unity.  In  this  limit  the  value  of  the  eigenvalues  k^  (R  =  1;  er  =  10)  a  are 
given  explicitly  by  the  quantities 


yfWr  k^(R  =0:er)</  =  l/VTo  k™(R  =0:  er  =  !())</. 


We  thus  observe  that  for  er  =  10  the  trajectories  start  for  R  «  1  at  the 
same  values  as  for  e,.  =  3,  but  decrease  to  values  for  R  ~  1  that  are  consid¬ 
erably  lower  than  for  er  =  3.  The  net  result  is  that,  compared  to  the  Er  =  3 
case,  for  er  =  10  the  trajectories  also  display  a  very  slow  decrease  for  in¬ 
creasing  R  initially.  This  behavior  persists  for  a  larger  range  of  R  as  n  in¬ 
creases.  A  clearly  defined  knee  is  immediately  evident  in  each  eigen¬ 
value  trajectory,  but  appears  earlier  with  R  for  er  =  10  compared  to  what 
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Figure  6.  TF.  eigenvalues  (R  :e,  =  10)  a  Tor  the  spherical  cavity  of  radius  a  Tilled  with  dielectric  e.  inter¬ 
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is  observed  in  figure  2  for  er  =  3.  The  £„  ,  (/?;  £r  =  10)  a  values  then  ex¬ 
hibit  a  rapid  decrease  over  a  much  larger  range  in  R  than  is  seen  in  figure 
2  for  er  =  3.  This  range  is  larger  the  larger  the  index  n  is.  We  note  in  fig¬ 
ure  6  that  this  rate  of  decrease  in  a  gradually  slows  as  R  increases. 
Finally,  as  R  increases  further,  the  rate  at  which  the  eigenvalues  decrease 
becomes  nearly  imperceptible.  This  latter  behavior  can  be  seen  to  persist 
for  a  shorter  range  of  R  as  the  index  n  increases.  Again  we  can  discern 
three  quite  distinct  regimes  in  the  eigenvalue  trajectories  of  figure  6, 
namely  those  of  a  small  inner  sphere,  an  intermediate  size  inner  sphere, 
and  a  large  size  inner  sphere. 

Results  obtained  from  equation  (32)  for  p  =  2  are  displayed  in  figure  7 
for  n  =  1  to  n  =  5.  Again  we  note  that  these  eigenvalue  trajectories  start 
at  R  =  0  exactly  where  the  corresponding  trajectories  for  er  =  3  start,  but 
they  end,  at  R  =  1,  at  values  V3 /VTo  times  those  for  er  =  3.  There  is  a 
considerable  amount  of  similarity  in  both  sets.  Examination  of  table 
2(b),  which  gives  numerical  values  for  the  k^iR'-Zr  =  10)  t/,  will  also 
show  the  behavior  of  these  eigenvalue  trajectories. 


Figure  7.  TE  eigenvalues  ^(R'Xr  =  10)  a  fur  the  spherical  cavity  of  radius  a  filled  with  dielectric  £,  inter¬ 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  10E,. 


Table  2a.  Numerical  Values  of  the  TE  Eigenvalues  k^(R\er  =  10)a  for  the  Spherical  Cavity 
Loaded  Internally  with  a  Dielectric  Sphere 


R«b/a 

n-1 

n=2 

n=3 

n«4 

n-5 

n=6 

0.05 

4.49328 

5.76346 

6.98793 

8.18256 

9.35581 

10.51284 

0.10 

4.48885 

5.76331 

6.98793 

8.18256 

9.35581 

10.51284 

0.15 

4.44690 

5.76035 

6.98773 

8.18255 

9.35581 

10.51284 

0.20 

4.19332 

5.72477 

6.98410 

8.18218 

9.35577 

10.51283 

0.22 

3.98261 

5.65198 

6.97555 

8.18118 

9.35565 

10.51282 

0.25 

3.63193 

5.30554 

6.87027 

8.16964 

9.35428 

10.51264 

0.27 

3.41190 

4.98023 

6.51613 

7.98705 

9.33367 

10.51108 

0.30 

3.11966 

4.52433 

5.90075 

7.22583 

8.50653 

9.75569 

0.40 

2.42757 

3.44379 

4.44930 

5.42871 

6.38311 

7.31778 

0.50 

2.00997 

2.79284 

3.57706 

4.35019 

5.10919 

5.85518 

0.60 

1.74458 

2.37268 

3.00662 

3.63843 

4.26405 

4.88225 

0.70 

1.57502 

2.09576 

2.62051 

3.14684 

3.67224 

4.19504 

0.80 

1.47448 

1.92305 

2.36881 

2.81429 

3.25991 

3.70540 

0.90 

1.42884 

1.83849 

2.23677 

2.62865 

3.01669 

3.40233 

0.95 

1.42199 

1.82475 

2.21361 

2.59360 

2.96743 

3.33677 

0.99 

1.42095 

1.82258 

2.20981 

2.58760 

2.95864 

3.32456 

Table  2b. 

Numerical  Values  of  the  TE  Eigenvalues  k 

%(*&  =  10)  a  for  the 

Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 

R=b/a 

n=1 

n-2 

n=3 

n*4 

n=5 

0.01 

7.72525 

9.09501 

10.41712 

11.70491 

12.96653 

0.05 

7.72328 

9.09499 

10.41712 

11.70491 

12.96653 

0.10 

7.60820 

9.09059 

10.41694 

11.70490 

12.96653 

0.12 

7.26468 

9.07331 

10.41605 

11.70485 

12.96653 

0.15 

6.28435 

8.76617 

10.40257 

11.70402 

12.96647 

0.17 

5.73488 

7.97978 

10.22983 

11.69839 

12.96609 

0.19 

5.35232 

7.22448 

9.32264 

11.38203 

12.96086 

0.20 

5.22407 

6.90371 

8.87423 

10.83897 

12.74741 

0.22 

5.07038 

6.39234 

8.09527 

9 . 86524 

11.60250 

0.25 

4.98980 

6.03912 

7.25222 

8.70091 

10.21430 

0.27 

4.97848 

5.98749 

7.09258 

8.24502 

9.47997 

0.30 

4.97672 

5.97476 

7.06823 

8.20943 

9.36395 

0.35 

4.90201 

5.94765 

7.05511 

8.20232 

9.36035 

0.37 

4.81017 

5.88405 

7.00951 

8.17039 

9.34000 

Table  2b.  (cont’d)  Numerical  Values  of  the  TE  Eigenvalues  k^  (/?;er  =  10)  a  for 
the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


Table  2c.  Numerical  Values  of  the  TE  Eigenvalues  k(^ 
ical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


(/?;£,  =  10)  a  for  the  Spher- 


R=b/a 


1 
5 
7 
0 
2 

4 

5 
.  1 7 
.19 
.20 
.22 
.25 
.27 
.30 
.35 
.39 


n**1 

n-2 

n-3 

n*4 

n-5 

10.90412 

12.32294 

13.69802 

15.03966 

16.35471 

10.89124 

12.32272 

13.69802 

15.03966 

16.35471 

10.79431 

12.32009 

13.69795 

15.03966 

16.35471 

9.44433 

12.24105 

13.69526 

15.03955 

16.35470 

8.52775 

11.29151 

13.66549 

15.03840 

16.35464 

8.24860 

9.93125 

12.65103 

15.01216 

16.35385 

8.21188 

9.52894 

1 1 . 84970 

14.44290 

16.35009 

8.19600 

9.31562 

10.66391 

12.77442 

15.01515 

8.19058 

9.29530 

10.49781 

11.75659 

13.44347 

8.17245 

9.29504 

10.49196 

11.73292 

12.98646 

8.04956 

9.28069 

10.49068 

11.72901 

12.97369 

7.55990 

9.04977 

10.42793 

11.71673 

12.97161 

7.15361 

8.63818 

10. 13765 

11.59456 

12.94049 

6.60346 

7.93799 

9.32125 

10.72313 

12.11663 

5.98979 

7.02022 

8. 13281 

9.27935 

10.43389 

5.81796 

6.65264 

7.57634 

8.54937 

9.55434 

5.80722 

6.61699 

7.51797 

8.48018 

9.49624 

5.78991 

6.58340 

7.46212 

8.42822 

9.46678 

5.70555 

6 . 53506 

7.42170 

8.37745 

9.38868 

5.59773 

6.43978 

7.32087 

8.24660 

9. 19462 

4.87694 

5.61702 

6.35145 

7.07965 

7.79857 

4.26651 

4.89482 

5.51222 

6.12089 

6.72171 

Table  2c.  (corn’d)  Numerical  Values  of  the  TE  Eigenvalues  (R;er  =  10)  a  for 
the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R«b/a 

n- 1 

n-2 

n**3 

n-4 

IA 

1 

C 

1 

0.80 

3.81758 

4.35949 

4.89019 

5.41234 

5.92759 

0.90 

3.52889 

4.00412 

4.46778 

4.92280 

5.37109 

0.95 

3.46157 

3.91555 

4.35652 

4.78772 

5.21127 

0.99 

3.44831 

3.89703 

4.33194 

4.75628 

5.17222 

Table  2d.  Numerical  Values  of  the  TE  Eigenvalues  ^  =  10)  a 

for  the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R*b/a 

n-1 

n«2 

n-3 

n«4 

0.01 

14.06618 

15.51460 

16.92362 

18.30126 

0.05 

14.00421 

15.51336 

16.92360 

18.30126 

0.10 

11.43270 

13.72226 

16.86368 

18.29972 

0. 12 

11.37161 

12.58012 

14.83143 

18.01644 

0.14 

11.36031 

12.52015 

13.78803 

15.53453 

0.15 

11.31392 

12.51957 

13.77198 

15.07830 

0.18 

10.59372 

12.36687 

13.75764 

15.06262 

0.20 

9.86784 

11.72783 

13.54654 

15.03454 

0.22 

9.32661 

10.90445 

12.69866 

14.54095 

0.25 

9.02345 

10.06308 

11.41212 

12.96304 

0.27 

9.00406 

9.89851 

10.96020 

12.18163 

0.29 

8.99128 

9.87796 

10.86449 

11.94701 

0.30 

8.95668 

9.87656 

10.85977 

11.93240 

0.34 

8.47016 

9.61876 

10.74045 

11.87308 

0.35 

8.29660 

9.45505 

10.60593 

11.77077 

0.40 

7.47526 

8.51649 

9.56963 

10.62991 

0.45 

6.90635 

7.76619 

8.66383 

9.57430 

0.49 

6.73440 

7.44741 

8.22620 

9.04498 

0.50 

6.72620 

7.41735 

8.17691 

8.98707 

0.55 

6.68379 

7.38512 

8.12963 

8.93470 

0.57 

6.59710 

7.32354 

8.07623 

8.87086 

0.60 

6.39671 

7.12858 

7.86922 

8.62393 

0.70 

5.63633 

6.28050 

6.91513 

7.54230 

0.80 

5.02010 

5.57962 

6.12794 

6.66742 
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Table  2d.  (cont’d)  Numerical  Values  of  the  TE  Eigenvalues  k{„ 4  ( R:er  =  10)  a 
for  the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R»b/a 


0.90 

0.95 

0.99 


n«*1 


4.59071 

4.47463 

4.44838 


n«2 


5.08176 

4.94037 

4.90650 


3 


5.56186 
5 . 39444 
5.35217 


n«4 


6.03332 

5.83934 

5.78794 


Table  2e.  Numerical  Values  of  the  TE  Eigenvalues 
k^  (R,tr  =  1 0)  ci  for  the  Spherical  Cavity  Loaded 
Internally  with  a  Dielectric  Sphere 


R=b/e 


R=b/a 


i 


tvc. 


For  er  =  10,  we  can  plainly  see  the  presence  of  a  second  knee  in  the 
curves.  Now,  however,  we  find  all  the  characteristics  are  strongly  accen¬ 
tuated.  We  see  in  figure  7  a  very  slow  decrease  in  k^l  (#;  Er  =  10)  a  as  R 
increases  from  zero.  This  behavior  persists  for  a  larger  range  in  R  as  n 
increases.  Now  the  first  knee  appears  for  Er  =  10,  earlier  than  for  Er  =  3. 
Thus,  for  example,  the  knee  appears  for  ^(R;  Er  =  10) a  at  R  =  0.2. 
whereas  for  k{J{(R\  Er  =  3)o,  it  is  at  R  =  0.3.  The  region  immediately 
following  the  first  knee  shows  a  much  more  rapid  decrease  as  R  increases 
for  Er  =  10  than  for  Er  =  3.  taking  place  over  a  smaller  range  of  R.  As  n 
increases,  this  behavior  is  accentuated  further.  As  we  can  see  in  figure  7. 
this  is  followed  by  a  region  of  quite  slow  decrease  in  the  eigenvalue  as  R 
increases  and  then  the  appearance  of  a  second  knee  which  is  sharper  for 
tr  =  10  than  for  Er  =  3.  Furthermore,  this  second  knee  occurs  at  lower  R 
values  for  the  larger  Er.  Figure  7  shows  that  the  second  knee  appears  at 
slightly  higher  values  of  R  for  larger  values  of  n.  After  rounding  the  sec¬ 
ond  knee  the  eigenvalue  decreases  rapidly  at  first  and  then  ever  more 
slowly  thereafter  as  R  increases.  The  larger  the  value  of  «,  the  more  pro¬ 
nounced  this  behavior  is.  Finally,  as  R  approaches  the  limit  unity  the 
trajectories  continue  to  decrease  as  R  increases,  but  very  slowly.  The 
range  of  this  latter  behavior  is  smaller,  but  the  rate  of  decrease  is  slower 
with  increasing  value  of  the  index  n. 

The  existence  of  three  regimes  is  again  quite  evident  in  figure  7.  For  the 
initial  range  of  R  for  which  the  eigenvalues  k ^  a  change  very  slowly  as  R 
increases,  we  have  a  small  interior  sphere.  For  the  upper  range  of  R  for 
which  the  eigenvalues  again  change  slowly  with  /?,  we  have  a  large 
dielectric  sphere  loading  the  cavity.  The  remainder  of  the  range  of  R 
from  the  onset  of  the  first  knee  then  corresponds  to  an  intermediate  inte¬ 
rior  dielectric  sphere.  For  the  intermediate  sphere,  as  R  increases,  a  vari¬ 
ety  of  behaviorial  aspects  are  displayed  by  the  system  eigenvalues. 

As  expected,  going  from  p  =  2  to  p  =  3  causes  a  third  knee  to  appear  in 
the  eigenvalue  trajectories  for  er  =  10,  just  as  we  saw  happen  for  e,.  =  3. 
Figure  8  shows  k^l (/? :  Er  =  10) «  for  n  =  I  to  n  =  5.  Table  2(c)  gives  nu¬ 
merical  values  for  these  eigenvalues.  Again  we  observe  a  sharpening  up 
of  the  characteristics  of  the  trajectories.  Starting  at  the  empty  cavity  reso¬ 
nances  for  R  =  0,  we  note  that  as  R  increases  the  k1^  (R:  Er  =  10)  a  values 
decrease  at  a  rate  slower  than  for  the  k[^(R.  Er  =  10)«  values,  but  for  a 


Figured.  I  K  eigenvalues  k'J}  (/?  f,  =  10)  a  fur  the  spherical  easily  of  radius  a  filled  with  dielectric  f  internally 
loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  IO0(. 


The  pattern  of  changes  induced  as  the  index  p  increases  by  unity  is  now 
evident.  Thus,  in  figure  10  we  show  only  the  n  =  \.  p  =  5  tra  jectory  w  ith 
an  additional  knee  and  shoulder  and  further  emphasis  of  the  detailed  be¬ 
havior  as  R  increases.  Table  2(e)  gives  numerical  values  for  the 
A,V(/?:  £r  =  10)  a  trajectory  shown  in  figure  9. 

It  should  be  quite  evident  what  one  would  expect  the  trajectories  for 
higher  order  eigenvalues  to  look  like  after  examining  the  given  sets  in 
detail.  It  is  a  relatively  simple  exercise  to  locate  with  good  accuracy  the 
limiting  values  for  R  =  1.  for  n  =  1  and  p  greater  than  5.  Thus,  for  p  =  5 
+  </.  say.  k\'\^  (R:er  =  10)  a  =  <5  +  q)  +  0.45.  This  result  can  be  deduced 
from  inspection  of  figures  6  through  10.  Similar  predictions  for  the  other 
values  of  n  can  be  made  in  the  same  wav.  which  means  we  can  find  these 
eigenvalues  (approximately)  without  directly  resorting  to  solving  equa¬ 
tion  (32). 

To  summarize,  some  very  intriguing  characteristics  are  exhibited  by  (he 
eigenvalue  trajectories.  Nevertheless,  we  shall  have  to  defer  the  explana¬ 
tion  of  these  to  a  subsequent  report  giving  the  detailed  behavior  of  the 
corresponding  spatial  distribution  of  the  field  components  in  the 
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Figure  10.  IF.  eigenvalues  k\'J(R:tr  =  10)  a  for  ihe  spherical  cavity  of  radius  a  filled  with  dielectric  £,  internally 
loaded  with  a  concentric  dielectric  sphere  of  radius  h  and  strength  E2  ~  ld£|- 

dielectric  loaded  cavity.  We  shall  have  to  be  content  here  with  merely 
noting  the  rich  variety  of  interesting  characteristics  displayed  by  the 
eigenvalue  trajectories. 

3.4  The  TE  Eigenvalues  for  er  =  1/3 

We  have  thus  far  examined  effects  on  the  eigenvalues  induced  by  loading 
the  interior  of  the  cavity  with  a  dielectric  sphere  whose  dielectric  constant 
is  larger  than  that  of  the  cavity,  i.e.,  er  =  e^/e,  >  1.  For  that  situation  we 
considered  representative  values  of  er  that  covered  a  substantial  range.  In 
an  attempt  at  completeness,  we  next  consider  effects  on  eigenvalues  of 
the  spherical  cavity  that  result  when  the  cavity  is  concentrically  loaded 
with  a  dielectric  sphere  whose  relative  dielectric  constant  is  less  than 
unity.  Many  of  the  essential  characteristics  can  be  seen  in  the  results  for 
Er  =  £-,/£,  =  1/3.  Therefore,  we  limit  ourselves  to  considering  only  this 
one  case.  Now,  when  R  =  0,  the  eigenvalues  y,a  are  once  again  the  same 
as  for  the  empty  spherical  cavity.  But,  in  the  other  limit,  as  R  approaches 
unity,  the  eigenvalues  are  multiplied,  rather  than  divided,  by  V3,  since 

\£  =  VT/3 . 


In  the  previous  two  cases  of  dielectric  loading,  we  solved  for  the  eigen¬ 
value,  y ,</.  knowing  that  we  would  be  able  to  find  the  corresponding 
wavenumber  in  the  denser  dielectric  region  by  means  of  the  equation 


If  we  wish  to  be  consistent  in  defining  our  eigenvalue  as  the  one  in  the 
region  of  smaller  dielectric  constant,  we  should  be  solving  for  y \a.  rather 
than  Y|</.  This  formulation  is  usually  desirable  when  the  medium  of 
lesser  dielectric  strength  is  air.  for  example.  But  y2a 's  equivalent  to  y^i 
multiplied  by  the  square  root  of  the  relative  dielectric  constant.  Please 
note  that  it  is  this  "normalized"  y^a  or.  equivalently,  ka.  that  is  listed  in 
the  tables  and  plotted  in  the  figures  which  relate  to  this  section.  The 
"true"  values  of  ka  are  achieved  by  multiplying  the  tabled  values  by 

VT/3 . 

Some  very  interesting  behavior  manifests  itself  immediately  in  figure  1 1. 
which  displays  the  eigenvalue  trajectories  k{„}  (R:  Er  =  1/3)  a  for  p  =  1 
and  n  =  1  to  n  =  6.  We  clearly  observe  that  there  is  a  considerable  range 
over  R  for  which  the  eigenvalues  increase  with  R .  but  at  an  astonishingly 
slow  rate.  This  behavior  clearly  becomes  even  more  pronounced  as  n  in¬ 
creases.  Thus,  the  total  change  for  n  =  1  as  R  goes  from  zero  to  0.3  is 
only  about  0.1.  whereas,  for  n  =  6  as  R  goes  from  0  to  0.55,  the  same  to¬ 
tal  change  is  incurred.  That  is.  for  n  =  6,  about  twice  the  range  in  R  re¬ 
sults  in  the  same  total  change.  Considered  from  the  point  of  view  of  the 
filling  factor  of  the  cavity,  it  is  quite  remarkable  to  note  that  substantial 
filling  of  the  cavity  by  the  dielectric  load  hardly  affects  the  eigenvalues 
k„l  a .  and  this  effect  persists  to  larger  R  as  n  increases.  Thus,  the  higher 
modes  n  >  I  appear  to  change  even  more  reluctantly  as  the  load  increases 
in  size. 

Following  this  range  of  slow  variation  the  eigenvalue  k\\ a  increases 
more  rapidly  with  R  at  first  and  then  exhibits  a  moderately  slow  rate  of 
increase  with  R.  Then  for  large  R  at  about  R  =  0.9.  the  eigenvalue  k\\}  a 
again  resumes  a  very  slow  rate  of  increase  as  R  increases  to  unity  where 
the  load  has  expanded  to  homogeneously  fill  the  cavity.  These  charac¬ 
teristics  are  more  enhanced  as  n  increases.  Indeed,  as  n  increases  a  larger 
range  of  variation  in  k^a  can  be  seen  in  the  intermediate  region  and  a 
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nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  E2  =  [/.?£,. 

shorter  range  in  the  final  slowly  varying  region  as  R  approaches  unity. 
Table  3(a)  contains  numerical  values  for  the  k^a  eigenvalues  of  figure 
11. 

Before  continuing  on  to  p  =  2.  let  us  compare  the  eigenvalue  trajectories 
for  p  =  I,  Er  =  3,  and  p  =  I .  £r  =  1/3;  i.e.,  the  curves  in  figures  2  and  1  I . 
Clearly  they  are  all  monotonic.  For  p  =  1  and  fixed  n  the  trajectories  for 
Er  =  1/3  display  increasing  values  of  the  eigenvalues  with  increasing  R. 
whereas  those  for  Er  =  3  sho»>  decreasing  values  as  R  increases.  Both 
sets  suggest  three  size  regimes  for  the  inner  dielectric  sphere — small,  in¬ 
termediate,  and  large.  The  small  inner  sphere  regime  appears  to  cover  a 
larger  range  of  sizes  for  Er  =  1/3  than  for  Er  =  3.  The  opposite  is  true  for 
the  large  inner  sphere  regime.  The  intermediate  regimes  are  about  the 
same.  Now  let  us  proceed  to  p  =  2  for  er  =  1/3. 

Figure  12  shows  the  p  =  2  eigenvalue  trajectories  for  n  =  1  to  n  =  5  and 
table  3(b)  contains  numerical  values  for  the  corresponding  eigenvalues 
(R :£r  =  1/3) a.  For  R  =  0,  they  start  out  where  expected,  then  show  a 
very  slow  increase  with  increasing  R  at  first.  This  behavior  persists  for  a 
smaller  range  of  R  than  we  saw  in  figure  11  forp  =  1  for  each  value  n. 


•».  % 


Table  3a.  Numerical  Values  of  the  TE  Eigenvalues 
Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


k{„\'  (Rxr  =  1/3)  a  for  the  Spherical 


R=b/a 


0.05 
0. 10 
0.20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
0.85 
0.90 
0.95 
0.99 


n=1 


2.59428 

2.59443 

2.59872 

2.62210 

2.68788 

2.82307 

3.06229 

3.45372 

4.00972 

4.26971 

4.42804 

4.48565 

4.49335 


n=2 


3.32753 

3.32754 
3.32813 
3.33555 
3.37083 
3.47099 
3.68799 
4.10280 
4.83103 
5.28293 
5.61867 
5.74678 
5.76333 


n=3 


4.03448 

4.03448 

4.03456 

4.03660 

4.05340 

4.12148 

4.30495 

4.70920 

5.52816 

6.14535 

6.71618 

6.95753 

6.98770 


n=4 


4.72420 

4.72420 

4.72421 
4.72473 
4.73207 
4.77548 
4.92383 
5.30154 
6.15939 
6.89590 
7.72475 
8.13244 
8.18219 


Table  3b.  Numerical  Values  of  the 
Spherical  Cavity  Loaded  Internally 


R=b/a  n=1  n=2 


TE  Eigenvalues  k^  (R'Xr  = 
with  a  Dielectric  Sphere 


5.40158 

5.40158 

5.40158 

5.40170 

5.40471 

5.43102 

5.54699 

5.89075 

6.75615 

7.57135 

8.64548 

9.27855 

9.35526 


1/3)  ci  for  the 


0. 

E'&S 

4.46025 

5.25101 

6.01433 

6.75783 

7.48623 

0. 

m 

4.46228 

5.25112 

6.01433 

6.75783 

7.48623 

0. 

20 

4.50698 

5.26123 

6.01619 

6.75813 

7.48627 

0 

30 

4.67880 

5.34694 

6.05175 

6.77099 

7.49046 

0 

40 

5.02509 

5.61047 

6.22883 

6.87885 

7.55093 

0 

50 

5.55227 

6.12289 

6.67413 

7.24173 

7.83216 

0 

60 

6.13356 

6.90351 

7.48865 

8.02226 

8.54929 

0 

65 

6.29287 

7.31436 

8.03667 

8.61560 

9. 14461 

0 

70 

6.32829 

7.56756 

8.56799 

9.32160 

9.92728 

0 

75 

6.33922 

7.61942 

8.83647 

9.91024 

10.78773 

0 

80 

6.45504 

7.64411 

8.86344 

10.07269 

11.23393 

0 

85 

6.81283 

7.88007 

8.97735 

10. 11151 

11.26856 

0 

90 

7.37109 

8.51928 

9.58584 

10.61248 

11.63328 

0 

95 

7.68393 

9.02563 

10.30954 

11.54703 

12.74418 

0 

99 

7.72494 

9.09451 

10.41636 

11.70383 

12.96506 

6.06959 

6.06959 

6.06959 

6.06962 

6.07079 

6.08608 

6.17424 

6.48127 

7.33493 

8.19897 

9.48357 

10.39934 

10.51206 


E«  I 


0.05  6.29589 
0.10  6.30579 
0.20  6.46361 
0.30  6.89422 
0.40  7.55005 
0.50  8.10967 
0.55  8.15397 
0.60  8.19431 
0.65  8.45088 
0.70  8.95862 
0.75  9.42703 
0.80  9.54680 
0.85  9.58831 
0.90  10.04004 
0.95  10.78049 
0.99  10.90325 


7.11466  7.90856 
7.11552  7.90861 
7.16979  7.92311 
7.46634  8.09436 
8.09855  8.62558 
8.99227  9.57439 
9.35481  10.15639 
9.45837  10.61068 
9.49536  10.71575 
9.82443  10.80233 
10.48337  11.35852 
10.89733  12.10373 
10.93208  12.25989 
11.23666  12.43592 
12.13873  13.43716 
12.32168  13.69629 


8.68315  9.44240 

8.68316  9.44240 
8.68637  9.44302 
8.77173  9.48058 
9.18115  9.77332 

10.08155  10.58411 
10.72741  11.22935 
11.43380  12.03077 
11.89237  12.83231 
11.94343  13.13739 
12.23209  13.22150 
13.13389  14.02286 
13.55397  14.79674 
13.64711  14.86768 
14.68489  15.88868 
15.03737  16.35175 


Table  3d.  Numerical  Values  of  the  TE  Eigenvalues  k(^  (R  ;er  =  1/3)  a 
for  the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R=b/a 

n*1 

n=2 

n=3 

n=4 

0.05 

8. 12245 

8.95740 

9.77086 

10.56624 

0.10 

6. 15202 

8.96111 

9.77119 

10.56626 

0.20 

8.49046 

9. 12083 

9.83003 

10.58388 

0.30 

9. 19020 

9.71212 

10.25741 

10.85562 

0.40 

9.90277 

10.66974 

11.19390 

11.68674 

0.45 

9.97521 

11.12443 

11.80829 

12.31948 

0.50 

10.04684 

11.28868 

12.36939 

13.05365 

0.55 

10.41778 

11.34393 

12.55668 

13.66935 

0.60 

10.98816 

11.77403 

12.65503 

13.79198 

0.65 

11.34156 

12.45693 

13.27072 

14.06253 

0.70 

11.37485 

12.75106 

13.98707 

14.95315 

0.75 

11.68947 

12.82383 

14.09741 

15.39889 

0.80 

12.46143 

13.52556 

14.51208 

15.56183 

0.85 

12.72179 

14.12917 

15.44163 

16.63609 

0.90 

12.84194 

14.19586 

15.54383 

16.87999 

0.95 

13.78510 

15.12991 

16.41936 

17.66522 

0.99 

14.06432 

15.51208 

16.92034 

18.29709 

At  the  upper  limit  of  /?,  a  short  range  of  R  is  shown  also  for  which  the 
eigenvalues  increase  very  slowly  with  R.  This  range  decreases  very  rap¬ 
idly  as  n  increases.  The  region  between  these  two  limits  reveals  the  pres¬ 
ence  of  an  indentation  towards  lower  ka  values.  The  range  in  R  covered 
by  this  "dent"  decreases  with  increasing  values  of  n.  It  also  produces  an 
intermediate  region  in  which  the  rate  of  change  of  (R  ;er  =  1/3)  a 
changes  from  fast  to  slow  to  fast  as  R  increases.  Although  the  region  of 
slow  variation  covers  a  smaller  range  in  R  as  n  increases,  it  simul¬ 
taneously  results  in  a  much  slower  rate  of  change  with  R. 

Compared  with  the  k^  a  trajectories  for  er  =  3  in  figure  3,  where  we  note 
that  the  eigenvalues  are  monotonic  and  decreasing  with  increasing  /?,  the 
eigenvalues  for  er  =  1/3  are  monotonic  and  increasing  with  increasing  R. 
The  curves  of  figure  12  for  er  =  1/3  do  not  display  the  sharp  characteris¬ 
tics  of  those  of  figure  3  for  er  =  3.  They  do,  however,  retain  many  of  the 
other  features  of  the  trajectories  of  figure  3. 

Figure  1 3  displays  the  eigenvalue  trajectories  for  p  =  3  and  for  n  =  1  to  n 
=  5.  Numerical  values  for  these  k{^  (R ;er  =  1/3)  a  curves  are  given  in 
table  3(c).  We  shall  first  consider  these  curves  with  regard  to  how  they 
vary  as  R  and  n  vary.  Next  we  compare  them  to  the  p  =  2  curves.  Then 
we  shall  make  some  comparisons  with  the  case  er  =  3,  i.e.,  for  the  situa¬ 
tion  in  which  the  dielectrics  in  the  cavity  are  interchanged.  Initially,  we 
observe  that  the  values  start  at  the  known  values  for  a  cavity  homoge¬ 
neously  filled  with  dielectric.  Since  we  have  the  stronger  dielectric  fill¬ 
ing  the  cavity,  these  R  =  0  values  for  each  n  are  reduced  by  the  factor 
VT ltr  =  V3,  as  we  discussed  earlier.  As  R  increases,  the  eigenvalue  for  n  = 
1  very  slowly  increases  with  R  up  to  about  R  =  0.1.  Further  increase  in  R 
causes  the  eigenvalue  to  increase  at  a  more  rapid  rate  until  about  R  =  0.5, 
where  again  it  dramatically  slows  its  rate  of  increase  until  R  has  increased 
to  about  R  =  0.6.  This  is  followed  again  by  a  more  rapid  increase  in 
k\ y  (R  ,er  =  1/3)  a  with  R  and  a  tapering  off  from  this  rate  of  increase  with 
further  increase  in  R,  starting  at  about  R  =  0.7.  We  observe  that  over  ap¬ 
proximately  the  range  0.5  5  R  5  0.7  there  is  present  in  the  k\]j  a  trajectory 
an  indentation. 


Immediately  beyond  this  range  in  R.  we  can  clearly  see  another  identa- 
tion  over  which  the  same  sort  of  behavior  in  k\\ a  is  clearly  visible.  Just 
after  this  second  indentation  in  the  trajectory  there  is  a  rapid  return,  as  R 
approaches  unity,  to  the  limiting  eigenvalue  for  the  cavity  homoge¬ 
neously  filled  with  dielectric. 

As  we  let  n  increase,  the  slowness  of  the  very  early  increase  in  k'^a  as  R 
increases  from  R  =  0  becomes  even  more  pronounced  and  persists  for  a 
larger  range  of  R.  Thus,  at  n  =  5  this  behavior  is  quite  obvious  in  figure 
13  up  to  about  R  =  0.3.  Clearly,  as  n  increases  further,  this  would  be 
even  more  enhanced  for  a  larger  range  in  R.  Also,  the  rapid  variation 
with  R  that  follows  this  initial  behavior  increases  further,  as  n  increases, 
over  a  slowly  increasing  range  of  R.  This  behavior  as  R  increases  is  then 
followed  by  a  short  range  of  R  over  which  k^a  increases  with  R,  but 
slowly.  As  n  increases,  this  region  of  the  trajectory  becomes  shorter  until 
at  n  =  5,  it  almost  is  undetectable  in  figure  1 3. 

Following  this  behavior,  we  find  that  k'^a  increases  more  rapidly  with  R 
through  an  indented  region.  We  note  in  figure  13  as  n  increases  this 
range  of  increase  with  R  grows  more  rapid.  Thereafter,  we  observe  the 
effecis  of  a  second  indentation  in  the  trajectories  and  a  repeat  in  the  be¬ 
havior  of  the  k'n'J  a  values  as  R  increases  through  the  range  corresponding 
to  the  identation  up  to  the  limit  at  R  =  1 . 

The  p  =  3  case  of  figure  1 3  primarily  differs  from  the  p  =  2  case  of  figure 
1 2  by  the  addition  of  one  more  identation  in  the  trajectories  and  the  resul¬ 
tant  characteristics.  If  we  compare  the  p  -  3  case  for  er  =  1/3  of  figure  1 3 
with  the  p  =  3  case  for  er  =  3  of  figure  8.  we  observe  the  former  is 
monotonically  increasing  with  R .  whereas  the  latter  is  monotonically 
decreasing  as  R  increases.  The  (R  \tr  =  1/3)  a  curves  are  considerably 
smoother  than  the  k^l(R:er  =  3 )ti  curves,  although  both  have  the  same 
number  of  indentations  in  the  same  sense  relative  to  ka.  Thus,  inter¬ 
changing  the  dielectrics  in  the  cavity  results  in  some  very  remarkable 
differences  in  the  corresponding  eigenvalue  trajectories. 

The  trends  in  the  TE  eigenvalue  trajectories  as  p  is  increased  further  are 
now  clearly  apparent.  We  see  in  figure  14  pretty  much  what  we  now  ex¬ 
pect  to  for  p  =  4.  This  is  also  evident  in  table  3(d)  for  the 
(Rxr  =  1/3)  a  numerical  values.  Clearly,  another  indentation,  the 
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third  one  now,  is  present  in  the  k^a  trajectories,  and  the  expected  behav¬ 
ior  is  exhibited  over  the  consecutive  ranges  of  R  as  we  go  from  R  =  0  to 
R=  I. 

We  have  just  completed  examining  a  large  number  of  sets  of  TE 
k„p  (R  :tr)  a  trajectories  for  each  value  of  p  and  for  a  substantial  range  of 
relative  dielectric  strength  er.  Next,  we  proceed  to  consideration  of  the 
corresponding  sets  of  TM  eigenvalues. 


I*  If  HSUOAlUfl  *2^1*  f,U 


Figure  14.  TK  eigenvalues  (/?  :tr  =  1/3)  a  fur  the  spherical  cavity  of  radius  a  filled  with  dielectric  E,  inter¬ 

nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  Ej  =  I/3E,. 


4.  The  Transverse  Magnetic  (TM)  Eigenvalues 


4.1  Introduction 

Just  as  we  did  in  section  3.1  earlier  for  the  TE  eigenvalues,  we  now  intro¬ 
duce  some  notation  conventions  for  the  TM  eigenvalues.  Using  the  iden¬ 
tities  just  preceding  equation  (32).  we  rewrite  equations  (29)  and  <301. 
respectively,  as  follows: 

£rln(y[trR  Yi  o  )  7V  ’  (a  ,Ro ) 

_  <  36 ) 

+  ['ItrRY\uin('JzrRy\a)] YJ'(Ra.a)  =0  . 


where 

(37) 


For  fixed  values  of  R  and  e,  we  consider  y,<;  the  variable  in  equation  (36) 
(just  as  in  the  TE  case).  The  values  Y,u  that  satisfy  this  equation  are  the 
eigenvalues  for  the  TM  modes.  Again  we  note  that  since  the  dielectrics 
w  ithin  the  cavity  are  lossless,  and  the  outer  cavity  wall  is  a  perfect  conduc¬ 
tor.  all  the  roots  of  equation  (36)  are  purely  imaginary.  Then,  as  for  the 
TE  eigenvalues,  for  each  value  of  the  index  n,  the  TM  eigenvalue  spec¬ 
trum  will  be  a  doubly  infinite  discrete  set  of  ordered  real  numbers,  each  of 
which  we  denote  by  k'^a  where  ik™a  is  a  root  of  equation  (36).  For 
Fixed  n  we  thus  have  k'^  a  >  k(^a  for  p'  >  p  where  p,p'  =  1 ,2,3... 

Each  value  of  kl£a  is  a  function  of  the  relative  dielectric  strength,  er,  and 
of  R,  the  ratio  of  the  radius  of  the  inner  dielectric  radius  6.  to  the  cavity  ra¬ 
dius  a.  Thus  k'^a  =  k'^(Rxr)a.  For  the  TM  case,  we  shall  use  the  same 
convention  for  the  dielectrics  as  in  the  TE  case  at  the  limits  R  =  0  and  R  = 
1,  or  equivalently,  b  -  0  and  h  =  a,  which  correspond,  respectively,  to  the 
absence  of  the  inner  dielectric  sphere  that  loads  the  cavity,  and  to  the  inner 
dielectric  load  sphere  completely  filling  the  cavity.  Thus,  we  can  immedi¬ 
ately  obtain  the  TM  eigenvalues  from  the  well-known  cavity  eigenvalues 
in  these  limits.  As  before,  for  Er  >  1.  the  eigenvalues  k~r  a  for  R  =  0  are 
the  usual  ones  for  the  spherical  cavity,  which  vary  with  R  until  R  =  I. 
whereupon  they  are  just  those  for  R  -  0  reduced  by  the  factor  \Cr.  For  •  r 
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<  1  they  start  at  R  =  0  as  the  reduced  normalized  values,  and  go  over  at  R  = 
1  to  the  unreduced  values.  The  eigenvalues  between  these  limits  must  be 
computed  using  equation  (36)  to  obtain  the  entire  picture.  This  has  been 
done,  and  we  now  proceed  to  consider  the  results  so  obtained  for  the  same 
range  of  dielectric  properties  as  in  the  TE  case. 

The  TM  Eigenvalues  for  er  =  3.00 

We  shall  first  consider  the  details  of  the  results  obtained  for  cr  =  3.  Here 
we  recall  that  the  inner  dielectric  sphere  of  radius  h  has  a  dielectric  con¬ 
stant  £,.  which  is  three  times  that  of  the  dielectric  constant  f ,  of  the  cavity 
dielectric  in  the  absence  of  the  inner  load. 

In  figure  15.  we  find  the  TM  eigenvalue  trajectories  k (R  x,  =  3)  u  for  p 
=  1  and  with  n  going  from  n  -  I  to  n  =  X  Table  4(a)  contains  numerical 
values  for  these  eigenvalues. 
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Table  4a.  Numerical  Values  of  the  TM  Eigenvalues  k'^(R:er  =  3)  a  for 
the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R-b/a 

n-1 

n-2 

n=3 

n«4 

n-5 

0.05 

2.74286 

3.87023 

4.97342 

6.06195 

7.14023 

0.10 

2.73706 

3.87007 

4.97342 

6.06195 

7.14023 

0.20 

2.69190 

3.86492 

4.97291 

6.06190 

7.14022 

0.30 

2.58539 

3.83249 

4.96545 

6.06035 

7.13991 

0.40 

2.43108 

3.73200 

4.92125 

6.04385 

7.13420 

0.50 

2.26204 

3.54208 

4.77545 

5.95374 

7.08509 

0.55 

2.17964 

3.42198 

4.64826 

5.84425 

7.00431 

0.60 

2.10058 

3.29330 

4.48936 

5.67699 

6.84888 

0.65 

2.02536 

3.16109 

4.30924 

5.45979 

6.60725 

0.70 

1.95399 

3.02859 

4. 11815 

5  21204 

6.30396 

0.60 

1.82162 

2.76782 

3.72715 

4.68780 

5.64199 

0.90 

1.69988 

2.50972 

3.32898 

4.15081 

4  97012 

0.95 

1.64157 

2.37657 

3.11535 

3.85667 

4.59870 

0.99 

1.59555 

2.26400 

2.92417 

3.58021 

4.23397 

Table  4a.  (cont’d)  Numerical  Values  of  the  TM 
Eigenvalues  k^(R:tr  =  3)<r  for  the  Spherical 

Cavity  Loaded  Internally  with  a  Dielectric  Sphere 

R-b/a 

n«6 

n=7 

n-8 

0.05 

8.21084 

9.27546 

10.33524 

0.10 

8.21084 

9.27546 

10.33524 

0.20 

8.21084 

9.27546 

10.33524 

0.30 

8.21078 

9.27545 

10.33524 

0.40 

8.20888 

9.27484 

10.33504 

0.50 

8.18396 

9.26270 

10.32929 

0.55 

8.13056 

9.22989 

10.31005 

0.60 

8.00051 

9.13009 

10.23846 

0.65 

7.74811 

8.88003 

10.00150 

0.70 

7.38917 

8.46338 

9.52239 

0.80 

6.58395 

7.50955 

8.41636 

0.90 

5.78326 

6.58770 

7.38191 

0.95 

5.33969 

6.07819 

6.81306 

0.99 

4.88643 

5.53809 

6.18923 

0.05 

6.10886 

7.44298 

8.72175 

9.96755 

11.18898 

0. 10 

6.05385 

7.43963 

8.72160 

9.96754 

11.18898 

0.20 

5.71699 

7.34454 

8.70488 

9.96504 

11 .  18864 

0.30 

5.33241 

6.92317 

8.49585 

9.89331 

11  .  16763 

0.35 

5.18844 

6.62768 

8.19692 

9.71259 

11.08859 

0.40 

5.05083 

6.34388 

7.79200 

9.31550 

10.81376 

0.45 

4.89193 

6.08469 

7.37505 

8.76706 

10.22011 

0.50 

4.70492 

5.83982 

7.00053 

8.22779 

9.51894 

0.55 

4.50438 

5.60029 

6.67709 

7.77279 

8.90497 

0.60 

4.30845 

5.36815 

6.39444 

7.40906 

8.42805 

0.65 

4.12893 

5. 15080 

6.14206 

7. 11422 

8.07612 

0.70 

3  97157 

4.95453 

5.91497 

6.86154 

7.80095 

0.80 

3.73125 

4.63748 

5.53332 

6.42664 

7.32360 

0.90 

3.59128 

4.42487 

5.24643 

6.06267 

6.87761 

0.95 

3.55456 

4.35349 

5.13403 

5.90342 

6.66576 

0.99 

3.53556 

4.30780 

5.05449 

5.78389 

6.50075 

0.05 

9.28897 

10.71235 

12.06358 

13.38012 

0. 10 

9.10628 

10.69302 

12.06224 

13.38005 

0.15 

8.77250 

10.57386 

12.04292 

13.37761 

0.20 

8.48286 

10.24119 

11.92761 

13.35179 

0.25 

8.24336 

9.77384 

11.54281 

13.19856 

0.30 

7.90934 

9.33457 

10.89296 

12.63995 

0.35 

7.49142 

8.90388 

10.27771 

11.761 78 

0.40 

7.13785 

8.47827 

9.78701 

11.07181 

0.45 

6.92502 

8.13346 

9.38185 

10.61430 

0.50 

6.83769 

7.91323 

9.06659 

10.25283 

0.55 

6.79459 

7.79473 

8.84867 

9.95403 

0.60 

6.69222 

7.69193 

8.68765 

9.71039 

0.70 

6.24923 

7.25848 

8.23221 

9. 18334 

Table  4c.  (cont’d)  Numerical  Values  of  the  TM  Eigenvalues 
k'nl  (R  ’Zr  =  3 )a  for  the  Spherical  Cavity  Loaded  Internally 
with  a  Dielectric  Sphere 


R-b/a 

n*1 

n-2 

n-3 

n«4 

0.80 

5.78179 

6.71995 

7.63368 

8.52910 

0.90 

5.47244 

6.33328 

7.17548 

8.00447 

0.95 

5.40210 

6.23206 

7.04158 

7.83653 

0.99 

5.38160 

6.19214 

6.97764 

7.74468 

Table  4d.  Numerical  Values  of  the  TM  Eigenvalues  k*‘, l  (R.er  =  3 )a  for 
the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 

R-b/a 

n»1 

n-2 

n-3 

n-4 

0.05 

12.41947 

13.91814 

15.31350 

16.67415 

0.10 

12.04265 

13.84993 

15.30689 

16.67365 

0. 15 

11.62708 

13.48050 

15.21499 

16.65837 

0.20 

11.24562 

12.84136 

14.72384 

16.49682 

0.25 

10.61792 

12.22503 

13.81468 

15.66128 

0.30 

10.15323 

11.59708 

13.07775 

14.54327 

0.35 

10.01586 

11.19032 

12.52571 

13.90196 

0.40 

9.91200 

11.02213 

12.18209 

13.44077 

0.45 

9.56698 

10.80230 

11.94798 

13.09572 

0.50 

9. 10805 

10.36373 

11.56787 

12.72224 

0.55 

8.75829 

9.90507 

1 1 . 06709 

12.21408 

0.60 

8.61046 

9.60926 

10.65945 

11.73258 

0.70 

8.45130 

9.41560 

10.34744 

11.26968 

0.80 

7.87959 

8.84116 

9.77566 

10.68876 

0.90 

7.37383 

8.25694 

9.12024 

9.96792 

0.95 

7.24289 

8.09214 

8.92173 

9.73616 

0.99 

7.21087 

8.04246 

8.85120 

9.64218 

V  V  *. 


!t  is  important  to  notice,  at  the  outset,  that  the  TM  eigenvalues  k\^ a  are  the 
lowest-lying  set  for  er  =  3.  Furthermore,  to  cover  the  same  frequency 
range  shown  in  figure  2  for  the  TE  values  k'^a.  we  require  TM  eigen- 
values  k~x  a  to  higher  order  (i.e..  to  n  =  7  compared  with  n  =  5  for  the  TE 
range).  This  circumstance  w  ill  be  seen  to  hold  true  for  higher  p  values. 

We  now  make  one  comparison,  i.e..  the  behavior  of  the  TM  trajectories 
relative  to  the  TE  trajectories  as  R  increases  from  zero  to  unity  or  as  the 
dielectric  load  increases  in  size  until  it  fills  the  cavity  completely.  Note 
then  that,  in  general,  the  TM  eigenvalue  trajectories  are  concave 
downward  for  the  entire  range  of  R,  whereas,  in  the  TE  case,  for  p  =  I.  the 
trajectories  change  over  at  about  R  =  0.5  from  concave  down  to  concave 
upward. 

The  TM  eigenvalue  trajectories  k\~j a  are  monotonic,  decreasing  as  R  in¬ 
creases.  Initially,  this  is  a  quite  slow  decrease  with  increasing  R.  for  a 
range  of/?  that  grows  with  increasing  value  of  the  order  index  n.  Indeed, 
we  can  observe  in  figure  15  that  for  n  =  1.  this  holds  up  to  about  /?=().! 
and  increases  so  that  at  n  -  X  it  holds  up  to  about  R  =  0.6.  Clearly,  this 
range  continues  to  grow  with  further  increase  in  n.  Following  this  initial 
behavior  the  eigenvalues  decrease  more  rapidly  with  increasing  R 
and  continue  in  this  manner  almost  to  R  =  1  itself.  Note  in  figure  15  that, 
as  n  increases,  the  eigenvalue  trajectory  shows  a  more  rapid  rate  of 
decrease  as  R  increases.  For  the  p  =  I  case,  we  cannot  clearly  discern  three 
distinct  regimes  corresponding  to  a  small-sized  inner  dielectric,  an  inter¬ 
mediate  one.  and  a  large  one. 

Before  continuing  to  higher  p  values,  we  pause  to  focus  on  a  very  remark¬ 
able  feature  exhibited  in  figure  15.  Notice  that  for  n  =  X.  we  find  that  k[]1  a 
changes  by  less  than  one  percent  (compare  table  4(a)),  as  R  increases  to  R 
=  0.6.  This  tells  us  that  we  have  an  inner  dielectric  sphere  whose  radius 
can  vary  from  zero  to  0.6  times  the  radius  of  the  entire  cavity,  with  very 
little  change  in  the  eigenvalue.  ky' j  a.  A  small  increase  in  the  size  of  the  in¬ 
ner  sphere  from  this  point  on.  then,  causes  a  rather  substantial  increase  in 
the  eigenvalue.  To  understand  this  interesting  property  requires  a  detailed 
study  of  the  corresponding  spatial  distributions  of  the  field  components. 
We  defer  this  discussion  to  a  subsequent  repo  t 
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for  (he  TM  p  =  2  and  TE  />  =  I  trajectories.  There  are  two  knees  present  in 
each  set.  The  TM  eigenvalues  k'~la  show  smaller  ranges  in  /?,  for  which 
we  have  a  very  slow  decrease  in  y^i  with  increasing  R,  than  we  found  for 
the  TE  case  of  k^ci.  The  TM  trajectories  then  reveal  a  succeeding  range 
over  which  a  decreases  more  rapidly  as  R  increases,  with  the  rate  in¬ 
creasing  with  //.  This  rate  of  variation  with  R  is  much  slower  than  we  ob¬ 
served  in  figure  3  for  the  TE  trajectories  for  k^a.  The  next  range  in  R 
shows  that  k{^ a  changes  slowly  with  increasing  R ,  but  not  nearly  as 
slowly  as  the  TE  trajectories  did  for  k^l  a.  For  the  remainder  of  the  range 
of  R,  both  k'fla  and  k'n'!  a  behave  about  the  same.  The  TM  trajectories, 
k{;l(RXr  =  3)t;,  still  can  be  thought  of  as  displaying  three  distinctly  sized 


regimes  for  the  size  of  the  loading  dielectric  sphere. 

Figure  18  gives  the  calculated  set  of  TM  eigenvalue  trajectories  for  p  =  4. 
Here  we  show  results  for  /;  =  1  to  n  =  4.  Numerical  values  of 
k'^(R  \£r  =  3)«  are  presented  in  table  4(d).  The  same  trend  we  saw  for  the 
TM  p  =  3  and  TE  p  =  2  eigenvalues  continues  here.  As  anticipated,  an 
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Figure  17.  I'M  eigenvalues  k(^(R\ fr  =  3)<j  fur  the  spherical  cavity  of  radius  a  filled  with  dielectric  6,  inter¬ 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  62  =  36,. 
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THE  TM  EIGENVALUES 

Figure  18.  TM  eigenvalues  k^}(R;er  =  3)  a  fur  the  spherical  cavity  of  radius  a  Tilled  with  dielectric  £,  inter¬ 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  3E,. 

additional  knee  appears  in  the  curves  for  the  k^a.  The  k^\a  and  a 
eigenvalue  trajectories  again  are  quite  similar,  with  the  TM  trajectories 
showing  much  less  pronounced  emphasis  of  the  features  over  the  succes¬ 
sive  ranges  of  R.  Of  course,  the  TM  eigenvalue  trajectories  for  k^a  are  at 
larger  ka  values  than  the  TE  trajectories  for  k1^  a. 

4.3  The  TM  Eigenvalues  for  er  =  10.00 

We  will  now  consider  the  TM  eigenvalues  for  the  inner  sphere  of  a  larger 
dielectric  constant  of  e2  =  l()e, .  First  we  examine  the  results  obtained 
from  equation  (36)  for  p  =  1  and  n  =  1  to  n  —  8.  These  results  are  the 
trajectories  plotted  in  figure  19  and  whose  numerical  values  of 
k%(R;er  =  1 0)  c?  are  given  in  table  5(a).  Again  we  begin  by  pointing  out 
that  for  all  /?,  the  TM  eigenvalues  k^  (R:Er  =  10)  a  lie  lower  than  the  TE 
eigenvalue  k\\  (R:er  =  10)  rr. 

In  fact,  the  TM  k\2' a  are  the  lowest-lying  eigenvalues  for  Er  =  10.  This  can 
be  seen  by  simply  comparing  the  curves  of  figures  6  and  19. 
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Figure  19.  TM  eigenvalues  k^(R\er  =  10)  a  for  the  spherical  cavity  of  radius  a  filled  with  dielectric  ,  inter¬ 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  10£,. 

We  have  a  rather  involved  situation,  in  which  there  are  a  fairly  large  num¬ 
ber  of  parameters  in  relation  to  which  the  different  eigenvalues  can  be 
compared;  i.e.,  TM  compared  with  TE,  er  =  10  relative  to  er  =  3,  and.  for 
given  er,  the  TM  eigenvalues  for  given  n  can  be  compared  for  different  p 
values,  or  for  different  n  values  all  other  parameters  kept  fixed.  In  spite  of 
this  complexity  we  shall  try  to  carry  through  the  comparisons  where  they 
appear  significant. 

Continuing,  we  compare  the  k^{R\zr  =  10)  a  trajectories  of  figure  19  with 
the  k„x  (R,£r  =  3) a  trajectories  of  figure  15.  For  n  =  1.  we  see  they  both 
begin  at  the  same  value  at  R  =  0.  As  R  increases,  they  both  decrease 
slowly  at  first  and  then  more  rapidly,  with  k\ j’  ( R  .zr  -  10)  a  having  its  rate 
of  change  somewhat  higher,  so  that  at  R  =  I,  we  have  the  ratio 
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Table  5a.  Numerical  Values  of  the  TM  Eigenvalues  k^(R:er  =  10)  a 
for  the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R=b/a 

n=1 

n=2 

n=3 

n=4 

0.05 

2.74214 

3.87023 

4.97342 

6.06195 

0. 10 

2.73105 

3.86993 

4.97341 

6.06195 

0.20 

2.64339 

3.86052 

4.97252 

6.06187 

0.30 

2.44378 

3.79690 

4.95825 

6.05896 

0.35 

2.31574 

3.71345 

4.92707 

6.04950 

0.375 

2.24897 

3.65122 

4.89390 

6.03630 

0.40 

2.18174 

3.57298 

4.83418 

5.99852 

0.45 

2.04872 

3.37092 

4.54925 

5.53510 

0.50 

1.92058 

3.13335 

4.14872 

4.99242 

0.60 

1.68455 

2.67539 

3.47315 

4.16142 

0.70 

1.47479 

2.29613 

2.96820 

3.55954 

0.80 

1.28336 

1.97686 

2.56689 

3.09291 

0.90 

1.09462 

1.67147 

2.19867 

2.68187 

0.925 

1.04438 

1.58649 

2.09550 

2.56951 

0.95 

0.99112 

1.49106 

1.97528 

2.43705 

0.975 

0.93313 

1.37652 

1.81883 

2.25449 

0.99 

0.89498 

1.29160 

1.68879 

2.08593 

Table  5a.  (cont’d  )  Numerical  Values  of  the  TM  Eigenvalues 
k^l  {R  .Er  =  10)  a  for  the  Spherical  Cavity  Loaded 

Internally  with  a  Dielectric  Sphere 

R=b/a 

n=5 

n=6 

n=7 

n=6 

0.05 

7.14023 

8.21084 

9.27546 

10.33524 

0.10 

7.14023 

8.21084 

9.27546 

10.33524 

0.20 

7.14022 

8.21084 

9.27546 

10.33524 

0.30 

7.13965 

8.21073 

9.27544 

10.33524 

0.35 

7.13700 

8.21002 

9.27526 

10.33519 

0.375 

7.  13233 

8.20844 

9.27472 

10.33501 

0.40 

7. 10384 

8. 13416 

9.06122 

9.96681 

0.45 

6.40689 

7.23980 

8.05514 

8.85948 

0.50 

5.76783 

6.51605 

7.24965 

7.97353 

0.60 

4.80582 

5.42951 

6.04110 

6.64448 

0.70 

4. 11466 

4.65122 

5.17663 

5.69448 

0.80 

3.58569 

4.06005 

4.52305 

4.97835 

0.90 

3.13486 

3.56829 

3.98876 

4.40023 
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R=b/a 

n=5 

n=6 

n-7 

n-8 

0.925 

3.01614 

3.44341 

3.85718 

4.26132 

0.95 

2.87747 

3.30045 

3.71001 

4.10939 

0.975 

2.68072 

3.09676 

3.50308 

3.90070 

0.99 

2.48199 

2.87605 

3.26741 

3.65560 

Table  5b.  Numerical  Values  of  the  TM  Eigenvalues  k'^  (/?:£,  =  10)  a  for  the 
Spherical  Cavity  Loaded  Internally  w  ith  a  Dielectric  Sphere 


R-b/a 

n«1 

n-2 

0.05 

6.10160 

7.44289 

0.10 

5.98960 

7.43678 

0.15 

5.69933 

7.39369 

0.20 

5.26552 

7.20470 

0.225 

5.02341 

6.94597 

0.23 

4.97386 

6.87218 

0.24 

4.87396 

6.70551 

0.25 

4.77340 

6.52213 

0.26 

4.67267 

6.33255 

0.30 

4.27969 

5.62072 

0.35 

3.84258 

4.92412 

0.39 

3.54785 

4.50996 

0.40 

3.48151 

4.42358 

0.41 

3.41790 

4.34339 

0.42 

3.35693 

4.26900 

0.45 

3.18869 

4.07635 

0.50 

2.95091 

3.82869 

0.53 

2.82963 

3.70623 

0.55 

2.75633 

3.63099 

0.60 

2.59564 

3.45729 

0.65 

2.46170 

3.29754 

0.70 

2.34909 

3. 14768 

0.80 

2.17202 

2.86919 

0.90 

2.03996 

2.60861 

0.95 

1.98493 

2.48136 

0.99 

1.94417 

2.37939 

n-3 

n«4 

n-5 

8.72175 

9.96755 

11.18898 

8.72148 

9.96754 

11.18898 

8.71724 

9.96717 

11 . 18896 

8.68543 

9.96260 

11 . 18834 

8.61675 

9.95149 

11.18654 

8.58611 

9.94672 

11.18578 

8.47633 

9.92863 

11. 18322 

8.27290 

9.85558 

11.17612 

8.01299 

9.58447 

11.06590 

7.00767 

8.33831 

9.61651 

6.04519 

7.15857 

8.24588 

5.47705 

6.44758 

7.41018 

5.36332 

6.30471 

7.24398 

5.26305 

6.18938 

7.16024 

5.17761 

6.11654 

7.14239 

5.00409 

6.03773 

7.12132 

4.83290 

5.94295 

7.06296 

4.73435 

5.85696 

6.98214 

4.66392 

5.78179 

6.88751 

4.47009 

5.52795 

6.48720 

4.25768 

5.21219 

6.02879 

4.03744 

4.88694 

5.60866 

3.60328 

4.28918 

4.89480 

3.18883 

3.75198 

4.28117 

2.98009 

3.47497 

3.95896 

2.80409 

3.22212 

3.63549 

Table  5e.  Numerical  Values  of  the  TM  Eigenvalues  k'~l  (R:er  =  10)  a  for  the 
Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R«b/e 

n«1 

n-2 

n-3 

n«4 

n-5 

0.01 

9.31621 

10.71301 

12.06359 

13.38012 

14.67012 

0.05 

9.26192 

10.71183 

12.06357 

13.38012 

14.67012 

0.10 

8.85483 

10.67402 

12.06114 

13.37999 

14.60711 

0.12 

8.55316 

10.60966 

12.05479 

13.37946 

14.67007 

0.14 

8.18958 

10.45346 

12.03635 

13.37748 

14.66988 

0. 15 

7.98515 

10.29326 

12.01612 

13.37513 

14.66961 

0.17 

7.54113 

9.67286 

11.87103 

13.36210 

14.66607 

0.18 

7.31408 

9.27880 

11.53410 

13.33712 

14.66588 

0.19 

7.09348 

8.89665 

11.01897 

13.09883 

14.65905 

0.20 

6.88492 

8.54797 

10.50791 

12.49429 

14.40867 

0.24 

6.21966 

7.62350 

8.97068 

10.45512 

12.02533 

0.25 

6.09741 

7.49808 

8.81040 

10.10740 

11.55074 

0.30 

5.69648 

7.06652 

8.51199 

9.88585 

11.16333 

0.32 

5.60586 

6.92199 

8.37712 

9.80510 

11.12933 

0.34 

5.53769 

6.78319 

8.20939 

9.66465 

11.04770 

0.35 

5.50855 

6.71509 

8.11276 

9.55969 

10.95706 

0.40 

5.35337 

6.36755 

7.53202 

8.73223 

9.87265 

0.45 

5.08304 

5.96893 

6.91103 

7.87507 

8.82188 

0.50 

4.73153 

5.54339 

6 . 35694 

7.17504 

7.98944 

0.55 

4.38900 

5.14612 

5 . 89083 

6.63380 

7.39156 

0.60 

4.08514 

4.79742 

5.50468 

6.23711 

7.07480 

0.65 

3.82306 

4.49806 

5.18430 

5.93650 

6.83130 

0.70 

3.59977 

4.24355 

4.91516 

5.67486 

6.55191 

0.80 

3.25674 

3.64864 

4.48193 

5. 18153 

5.91635 

0.90 

3.04264 

3.57702 

4.13069 

4.70382 

5.28066 

0.95 

2.98287 

3.47554 

3.97044 

4.46826 

4.96624 

0.97 

2.96632 

3.43930 

3.90782 

4.37386 

4.83778 

0.99 

2.95251 

3.40466 

3.84574 

4.27893 

4.70625 

Table  5d.  Numerical  Values  of  the  T\1  Eigenvalues  (/?:e,  =  10)  a 
for  the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 


R»b/a 

n-1 

n-2 

fl«3 

n«4 

0.01 

12.48495 

13.92052 

15.31356 

16.67415 

0.05 

12.34864 

13.91620 

15.31346 

16.67415 

0  09 

11.66566 

13.83091 

15.30744 

16.67380 

0.10 

11.39210 

13.75563 

15.30059 

16.67324 

0.11 

11.07169 

13.61101 

15.28721 

16.67202 

0.12 

10.70612 

13.30020 

15.25973 

16.66941 

0.  13 

10.32192 

12.73787 

15.18674 

16.66379 

0.14 

9.95744 

12.08968 

14.81769 

16.64928 

0.15 

9.64061 

11.52244 

13.99140 

16.51772 

0.16 

9.38302 

11.10856 

13.18860 

15.62409 

0.18 

9.03329 

10.66961 

12.17109 

13.93284 

0.20 

8.83911 

10.40660 

11.97936 

13.37331 

0.22 

8.72520 

10.17569 

11.82085 

13.31414 

0.25 

8.56776 

9.83956 

11.43985 

13.10634 

0.27 

8.37583 

9.58818 

11.05765 

12.71313 

0.28 

8.23946 

9.44241 

10.83503 

12.41234 

0.30 

7.91356 

9.10830 

10.36626 

11.75247 

0.35 

7.09866 

8.22480 

9.32598 

10.41587 

0.40 

6.54310 

7.56107 

8.64055 

9.79060 

0.45 

6.32088 

7.20096 

8.21410 

9.33433 

0.50 

6.23596 

7.00887 

7.88291 

8.83148 

0.55 

5.99439 

6.73500 

7.49604 

8.28084 

0.60 

5.63298 

6.34331 

7.04389 

7.74185 

0.70 

4.95085 

5.58490 

6.20868 

6.83158 

0.80 

4.42896 

4.99928 

5.57003 

6.15672 

0.90 

4.07959 

4.59918 

5. 12427 

5.66306 

0.95 

3.98594 

4.47725 

4.96790 

5.46056 

0.99 

3.95318 

4.41489 

4.86623 

5.30972 
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Thus,  all  the  k1^  (Rxr  =  10)  a  and  (Rxr  =  3 )a  will  have  the  same  ratio 
at  R  =  1  for  each  n. 

For  n  =  1,  both  trajectories  look  pretty  much  alike  for  er  =  10  and  er  =  3. 
However,  we  have  already  noted  that  as  R  increases  the  kn\  a  values 
decrease  faster  for  er  =  10  than  for  er  =  3.  There  are  some  other  dis¬ 
similarities  that  arise  when  we  compare  figure  l1)  with  figure  15.  Clearly, 
for  er  =  10.  the  trajectories  decrease  slowly  with  increasing  R\  the  rate 
decreases  as  «  increases.  This  was  also  true  of  the  eigenvalue  trajectories 
k'*{  ( RXr  =  3 )a.  However,  at  the  knee  we  observe  a  much  sharper  turning 
for  the  er  =  10  trajectories  relative  to  those  for  er  =  3.  One  new  feature  is 
that  as  n  increases,  for  £r  =  10.  we  have  a  decreasing  initial  range  of  R  over 
which  this  very  slow  changing  of  k{^a  occurs.  In  addition  to  this  be¬ 
havior.  we  observe  that  the  k'^(Rx,  =  1 0) trajectories  in  figure  19,  just 
beyond  the  knee,  are  clearly  concave  upward.  This  behavior  persists  for  a 
large  range  of  R  and  is  immediately  followed  by  a  region  of  rapid  decrease 
in  a,  as  R  goes  to  unity.  In  figure  15,  there  is  only  a  hint  of  this  con¬ 
cavity  in  the  k'~l  a  trajectories  for  er  =  3.  Note  that  the  k(*j  (Rxr  =  10)  a 
trajectories  somewhat  resemble  the  k[\(RXr  =  10)  u  trajectories  (compare 
fig.  6  and  19),  except  near  R  =  I . 

In  summary,  then,  our  considerations  thus  far  have  been  to  compare  the 
TM  k'^(RXr  =  10)  a  to  the  TE  k'^(RXr  =  10)  a  and  also  to  compare  the 
TM  k'£(RXr  =  10)  to  the  TM  k**{RXr  =  3)<j. 

Next,  we  examine  the  characteristics  of  the  TM  eigenvalue  trajectories  for 
p  =  2,  k^  {RXr  =  10)  (i.  when  n  goes  from  1  to  5.  The  k„la  curves  are 
shown  in  figure  20,  and  numerical  values  are  given  in  table  5(b).  Compar¬ 
ing  figure  20  with  figure  7,  we  see  that  there  are  two  knees  in  all  the  TE 
trajectories  k(Jj  a,  whereas  only  for  n  £  4  do  we  see  two  clearly  defined 
knees  in  the  TM  eigenvalue  curves  for  k'^ ( RXr  =  10) a.  For  n  >  4.  the 
TE  and  TM  trajectories  are  similar,  except  for  values  of  large  R.  For  n  < 
3.  the  k[^(RXr  =  10)  a  trajectories  resemble  more  closely  the  TE  eigen¬ 
value  trajectories  k1^  (Rxr  =  10)  «.  There  now  appears  to  be  a  transition 
of  some  sort  occurring  in  the  inflection  characteristics  of  the  curves  as  we 
vary  n 

Going  to  higher  order  in  p.  we  can  follow  the  development  of  this  transi¬ 
tion.  Thus,  we  consider  next  the  p  =  3  TM  eigenvalues  k^u.  The 


CONDUCTING  SPHERE  OE  RADIUS  a 


Figure  20.  TM  eigenvalues  k(^  (R'.fr  =  10)  a  for  the  spherical  cavilv  i»f  radius  a  filled  with  dielectric  E1  inter¬ 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  E,  =  I  Of,. 

trajectories  are  shown  in  figure  21  for  n  =  1  to  n  -  5.  and  in  table  5(0  are 
corresponding  numerical  values.  We  first  compare  the  TE  trajectories 
(R.Sr  =  10)  o  and  the  TM  trajectories  k^(R;£r  =  10)  a.  All  the  TE 
trajectories,  as  can  be  seen  in  figure  X.  have  increased  the  number  of  knees 
to  three  for  all  values  of  n.  Figure  21  reveals  two  well-defined  knees  oc- 
curing  in  the  TM  trajectories  for  values  of  n  <  4.  The  n  =  5  TM  trajectory, 
k^x(R:e,  =  10)</,  shows  just  a  vestige  of  a  third  knee,  but  at  a  rather  large 
value  of  R. 

Perhaps  as  n  increases  much  further,  the  k'^a  will  display  a  clearly  de¬ 
fined  third  knee.  If  we  compare  the  k'^a  trajectories  with  the  TE  trajec¬ 
tories  k\ V  (RXr  =  10)  a  of  figure  7.  we  observe  that  for  n  =  1,2,  and  3  the 
TM  trajectories  appear  to  more  or  less  resemble  the  corresponding  TE 
trajectories.  Considerable  differences  are  obvious.  Again,  we  appear  to 
find  that  a  transition  in  the  TM  trajectory  characteristics  is  taking  place. 

Comparing  figure  17  for  the  k^(R:er  =  3 ) </  with  figure  21  for  the 
k™(RXr  -  10)  we  see  that  for  n  =  I  and  n  =  2.  the  trajectories  are 
similar.  For  the  larger  values  of  n  this  is  no  longer  the  case. 
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THE  TM  EIGENVALUES  fT)a 

Figure  21.  TM  eigenvalues  A („:,’(/?;er  =  1 0)  </  fur  the  spherical  cavity  of  radius  a  filled  with  dielectric  E,  inter* 
nallv  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  IOE,. 

For  p  =  4.  we  find  that  the  same  trends  continue  for  the  TM  trajectories 
k'~l  (R  .Er  =  10)  a.  Figure  22  shows  these  trajectories  for  n  from  n  -  I  to  n 
=  5.  Table  5(d)  contains  corresponding  numerical  values  for  the  eigen¬ 
values  k'„\ '(/?:£,  =  10) a.  Let  us  start  by  first  making  the  observation  from 
figure  22  that,  for  R  from  about  0.5  up  to  unity,  all  the  TM  trajectories 
manifest  the  more  familiar  behavior  as  R  increases  that  we  had  become  ac¬ 
customed  to.  Furthermore,  there  are  three  distinct  knees  in  each  curve.  In 
fact,  all  the  familiar  features  are  present  in  each  trajectory.  At  this  value  of 
p,  it  seems  that  the  transitioning  in  the  behavior  of  the  eigenvalue  trajec¬ 
tories  has  pretty  much  been  completed.  Once  again  we  notice  upon  com¬ 
paring  figure  22  with  figure  8  that  for  the  same  er  the  TM  eigenvalues  for 
the  index  p  are  very  similar  to  the  TE  eigenvalues  belonging  to  the  index 
one  lower  in  value. 

Although  we  have  not  calculated  k^(R:zr  =  10)  a  for  p  >  4.  we  anticipate 
that  the  trend  will  continue,  and  the  familiar  pattern  of  behavior  will  result. 

We  have  now  seen  TM  eigenvalues  for  a  substantial  range  of  £r  >  1  and 
the  corresponding  effects  due  to  the  strength  of  the  dielectric  load. 
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THE  TM  EIGENVALUES  h^'(R  (f)t 

Figure  22.  TM  eigenvalues  (/?;fr  =  IO)<i  for  the  spherical  cavity  of  radius  a  filled  with  dielectric  E,  inter¬ 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  h  and  strength  E2  =  IOE,. 

4.4  The  TM  Eigenvalues  for  Er  =  1/3 

The  final  TM  eigenvalue  trajectories  we  shall  consider  are  those  for  which 
we  interchange  the  dielectric  strengths  within  the  cavity.  We  now  have  the 
case  of  Er  =  1/3.  Figure  23  shows  the  first  set  of  these  eigenvalue  trajec¬ 
tories  corresponding  to  p  =  I.  namely,  the  k'^  ( R:er  =  1/3)  a  for  n  =  1  to  n 
=  9.  Table  6(a)  contains  corresponding  numerical  values. 

Note  first  that  for  all  values  of  R  the  TM  eigenvalues  k{^(R:zr  =  1/3)  a  lie 
lower  than  the  lowest  TE  eigenvalues  for  Ef  =  1/3.  which  are  the 
kW'iRie,  =  1/3)  a. 

Just  as  for  the  TE  trajectories  =  1/3) a  (compare  fig.  1 1 ).  the  TM 

trajectories  of  figure  23.  k'^(R :er  =  1/3)  </.  display  the  very  remarkable 
feature  of  increasing  very  slowly  at  first  as  R  increases  from  R  =  0,  and 
changing  even  more  slowly  w  ith  R  as  n  increases,  for  an  increasing  range 
of  R.  Thus,  for  example,  for  n  =  9.  k'^'a  changes  by  considerably  less  than 
one  percent  as  R  changes  from  R  =  0  to  about  R  =  0.65.  This  is  indeed  a 
very  noteworthy  feature:  the  radius  of  the  inner  sphere  increases  to  about 
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R»b/a 


0.05 
0.  10 
0.20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 


n-1 


1.58442 

1.58674 

1.60414 

1.64733 

1.72508 

1.84634 

2.01985 

2.24474 

2.48256 


2 


2.23449 

2.23456 

2.23682 

2.24998 

2.28995 

2.37795 

2.54163 

2.81645 

3.21933 


3 


2.87141 

2.87141 

2.87165 

2.87488 

2.89205 

2.94671 

3.07754 

3.34626 

3.84044 


n=4 


3.49987 
3.49987 
3.49989 
3 . 50060 
3.50723 
3.53842 
3 . 63556 
3.87599 
4.40321 


n«5 


4.12241 

4.12241 

4.12241 

4.12256 

4.12494 

4.14174 

4.21074 

4.41647 

4.94112 
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^  Table  6a.  (corn’d)  Numerical  Values  of  the  TM  Eigenvalues  k™  (/?:£,  =  1/3)  a 

$ 

& 

||  for  the  Spherical  Cavity  Loaded  Internally  with  a  Dielectric  Sphere 

» 

'S'j 

fi-b/a 

n»1 

n«2 

n-3 

n«4 

n-5 

X' 

j£  0.85 

2.58072 

3.43959 

4.17773 

4.82901 

5.42371 

W  0.90 

2.65484 

3.63064 

4.51880 

5.32948 

6.07207 

Si 

ja  0.95 

2.70661 

3.77138 

4.78887 

5.76744 

6.71068 

W  0.99 

iS 

2.73699 

3.85265 

4.94139 

6.01218 

7.06961 

K 

Table  6a.  (corn  d)  Numerical  Values  of  the  TM  Eigenvalues 

•V 

=  1/3)  a  for  the  Spherical  Cavity  Loaded 

r'j 
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R-b/a 

n*1 

n«2 

n«3 

n«4 

n»5 

0.60 

4.61418 

5.58355 

6.30254 

6.88120 

7.41095 

0.70 

4.74989 

5.83512 

6.83170 

7.67987 

8.38242 

0.75 

4.91025 

5.96067 

7.00539 

7.98584 

8.86366 

0.80 

5.18611 

6.18401 

7.20685 

8.22912 

9.22090 

0.85 

5.57252 

6.58177 

7.57069 

8.56728 

9.57072 

0.90 

5.92529 

7.08509 

6.16141 

9. 18788 

10.18974 

0.95 

6.08487 

7.37546 

8.60637 

9.79266 

10.94311 

0.99 

6.11428 

7.43672 

8.71029 
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11 . 16424 

Table  6c. 

Numerical  Values  of  the  TM  Eigenvalues 

=  1/3)  a  for  the  Spherical  Cavity  Loaded 

Internally  with  a  Dielectric  Sphere 

R=b/a 

n=6 

n=7 

n=8 

0.05 

7.15422 

7.83968 

8.51786 

0.10 

7.15422 

7.83968 

8.51786 

0.20 

7.15424 

7.83968 

8.51786 

0.30 

7.15647 

7.84029 

8.51802 

0.40 

7.19311 

7.85784 

8.52595 

0.50 

7.38866 

7.99383 

8.61649 

0.60 

7.93396 

8.46664 

9.01445 

0.70 

8.98304 

9.52724 

10.04568 

0.75 

9.62647 

10.28862 

10.87838 

0.80 

10.15707 

11.01921 

11.79863 

0.85 

10.57054 

11.55457 

12.51114 

0.90 

11.18270 

12.17409 

13.16565 

0.95 

12.06366 

13.15887 

14.23260 

0.99 

12.35869 

13.53702 

14.70191 

Table  6c. 

(cont'd)  Numerical  Values  of  the  TM  Eigenvalues 

k'£(**r 

=  1/3)  a  for  the  Spherical  Cavity  Loaded  Internally 

with  a  Dielectric  Sphere 

R«b/a 

n=l 

n-2 

n*=3 

n=4 

0.05 

6.38931 

6.18546 

6.96492 

7.72502 

0.10 

5.44748 

6.19321 

6.96554 

7.72506 

0.20 

5.72571 

6.32290 

7.00705 

7.73533 
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Table  6c.  (coin'd)  Numerical  Values  of  the  TM  Eigenvalues 
krf  ( R :er  =  1/3)  a  for  the  Spherical  Cavity  Loaded  Internally 
with  a  Dielectric  Sphere 


R=b/a 

n=1 

n-2 

n-3 

n*4 

0.30 

6.11831 
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7.24937 
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7.93460 
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11.10187 

0.80 
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0.85 

8.14174 

9.37617 

10.61980 

11.85040 

0.90 

8.71535 

9.87300 

11.00993 

12.14933 

0.95 

9.23482 

10.57775 

11.86018 

13.09449 

0.99 

9.31459 

10.70821 

12.05511 

13.36718 

Table  6c.  (cont’d)  Numerical  Values  of  the  TM 

Eigenvalues  k[^  (/?;£,  = 

1/3)  a  for  the  Spherical 

Cavity  Loaded  Internally  with  a  Dielectric  Sphere 
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n=5 

n=6 
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0.30 
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9.22811 
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0.40 
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0.45 
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0.50 
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0.55 
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11 . 15398 

0.60 
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0.65 
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0.70 
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0.75 
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0.80 
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0.90 
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0.10 
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8.20645 
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10.02946 
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10.55160 
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9.10309 

9.86534 

10.78236 
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12.41425 

12.99682 
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9.52710 

10.31643 

11.09718 

11.99178 

12.88249 

13.64077 

9.74453 

10.83041 

11.63100 

12.40714 
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14.15221 

0.65 

9.78081 

11.09210 

12.20129 

13.06583 

13.84316 

14.65510 

0.70 

10.00693 

11.16375 

12.41521 

13.59417 

14.60022 

15.45226 

0.75 

10.61125 

11.56821 

12.60217 

13.74771 

14.93253 

16.07729 

0.80 

11.08410 

12.30813 

13.35163 

14.31296 

15.29773 

16.34951 

0.85 

11.14081 

12.53845 

13.87439 
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16.29569 

17.36492 

0.90 

11.41012 

12.69611 

13.98595 

15.27241 

16.54870 

17.80976 

0.95 

12.29011 

13.63758 

14.92764 

16.17285 

17.38395 

18.57040 

0.99 

12.48354 

13.91571 

15.30558 

16.66232 

17.99221 

19.29977 
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2/3  of  the  cavity  radius  with  almost  no  change  in  the  eigenvalue  of  the  sys¬ 
tem.  As  is  to  be  expected,  as  R  increases  beyond  this  range,  k^a  in¬ 
creases  more  rapidly  with  R.  This  rate  of  increase  is  amplified  further  with 
increasing  n.  Although  this  is  quite  similar  to  the  behavior  of  the  TE 
eigenvalues  k^}(Rxr  =  l/3)</  as  R  becomes  nearly  unity,  the  two  sets  of 
trajectories  display  deviant  behavior.  Whereas  the  ti  slow  up  abruptly 
for  the  upper  range  of  /?,  the  k'*{  a  appear  to  continue  smoothly  and 
gradually. 

Increasing  p  to  2,  we  obtain  from  equation  (36)  the  TM  eigenvalues 
k^l(Rxr  =  1/3) ti,  whose  trajectories  for  each  value  of  the  index  n  are 
shown  in  figure  24,  where  n  goes  from  n  =  1  to  n  =  9.  Corresponding  nu¬ 
merical  values  are  given  in  table  6(b).  Comparison  of  the  TM  trajectories 
of  figure  24  with  those  of  the  k'~j  a  in  figure  23  reveals  that  for  p  =  2  a  dent 
appears  in  each  curve.  The  dent  appears  in  the  midrange  of  R  for  n  =  1. 
As  n  increases,  the  range  of  R  decreases  over  which  the  dent  is  present. 
Furthermore,  the  dent  moves  to  higher  R  values:  that  is.  it  is  present  for 
larger  inner  spheres,  as  n  increases.  Finally,  we  observe  in  figure  24,  at 
the  higher  values  of  n,  the  dent  in  the  trajectories  k^l  a  seems  to  be  wash¬ 
ing  out.  Again  comparing  the  eigenvalue  trajectories  in  figure  24  with 
those  in  figure  23  we  see  that  prior  to  the  region  of  the  dent,  i.e..  for  lower 
R  values,  the  curves  are  quite  similar.  Flowever,  the  initial  range  of  R,  for 
which  the  k^a  values  seem  to  change  only  very,  very  slowly  with  increas¬ 
ing  R.  decreases  substantially  compared  to  the  corresponding  range  for  the 
k^l  a  trajectories.  At  the  upper  end  of  the  range  of  R.  we  note  that  the 
*£’(/?;£,  =  1/3) «  trajectories  increase  slightly  faster  than  the 
k{„ (R  xr  =  1/3)  a  curves. 

Consider  next  the  p  =  3  TM  eigenvalue  trajectories,  (Rxr  =  1/3)  a,  for 
n  =  1  to  n  =  7  shown  in  figure  23.  Table  6(c)  gives  numerical  values  for 
the  eigenvalues  k$(Rxr  -  1/3) a.  A  second  indentation  is  now  seen  in 
each  trajectory.  For  n  =  1,  the  dents  appear  at  about  one-third  and  two- 
thirds  of  the  range  of  R.  As  n  increases,  these  dents  move  to  higher  R 
value,  with  the  lower-lying  one  washing  out  rather  rapidly  as  n  increases. 
We  also  observe  that  the  initial  range  of  /?,  over  which  each  trajectory  dis¬ 
plays  a  very  slow  increase  with  increasing  R,  shortens  up  as  p  increases.  It 
is  also  of  interest  to  compare  the  trajectories  for  ( Rxr  =  1/3)  a  in  figure 
23  with  the  TE  trajectories  for  k^(Rxr  =  1/3)  a  in  figure  13.  Both  sets 
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have  two  indentations.  Indeed,  both  sets  are  very  similar  in  detail.  Note 
that  the  TM  eigenvalue  trajectories  tend  to  lose  the  sharpness  of  the  inden¬ 
tation  characteristic,  compared  to  those  in  the  TE  trajectories  of  figure  1 3. 


The  final  set  of  eigenvalue  trajectories  calculated  from  equation  (36)  that 
we  consider  is  the  curves  plotted  in  figure  26  for  p  =  4.  i.e.,  the  TM  trajec¬ 
tories  k{^(R:er  =  1/3) a  for  n  =  1  to  n  =  6.  Numerical  values  for  these 
eigenvalues  are  in  table  6(d).  Each  eigenvalue  trajectory  now  has  three  in¬ 
dentations,  each  successively  shifted  upward  in  R  in  quarterly  intervals  in 
R.  Just  as  we  saw  for  the  lower  values  of  the  index  p ,  the  distinctness  of 
these  indentations  in  the  trajectories  fades  out  as  n  increases,  with  the 
lower  ones  fading  faster  for  a  given  value  of  n.  Again  we  observe  that  the 
initial  slow  increase  with  R  for  the  early  range  of  R  has  decreased  as  we 
increase  p.  We  further  see  that  increasing  the  index  p  results  in  predictable 
changes  in  the  trajectories.  Finally,  we  can  compare  the 

k%(Rxr  =  1/3  )a  trajectories  of  figure  26  with  the  TE  k'^(R:er  =  1/3)  a 
trajectories  of  figure  14.  Clearly  they  look  very  much  alike.  Notice  that 
the  TE  trajectories  are  shifted  up  in  ka  relative  to  the  TM  trajectories. 


Figure  24.  TM  eigenvalues  k^iR'-Cr  =  1/3) <i  for  the  spherical  cavity  of  radius  a  filled  with  dielectric  E,  inter¬ 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  h  and  strength  E2  =  1/3E,. 
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Figure  25.  TM  eigenvalues  ljl2,)(/?;er  =  l/3)u  for  the  spherical  cavity  of  radius  a  filled  with  dielectric  £,  inter 
nally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  E2  =  I/3E,. 


CONCENTRIC  DIELECTRIC  SPHERE 
RADIUS  b 

DIELECTRIC  CONSTANT  lt  -  111c, 


Further  increasing  the  index  will  add  additional  indentations  to  the  trajec¬ 
tories.  Otherwise,  the  trends  already  discerned  will  merely  continue  to 
manifest  themselves.  This  will  he  the  case  if  the  dielectric  system  is  the 
same.  Changing  the  relative  dielectric  strengths  w  ith  er  <  I  may  introduce 
some  new  variations  in  the  eigenvalue  trajectories.  Since  we  have  not  in¬ 
vestigated  this  aspect  of  the  eigenvalue  problem  for  the  generic  system 
studies,  we  shall  not  go  any  further  into  this. 


and  Discussion 


General  Summary 

We  have  examined  a  very  large  number  of  eigenvalues  calculated  from 
equations  (32)  and  (36).  We  shall  now  summarize  briefly  the  general 
characteristics  we  have  found  for  these  eigenvalues.  First,  for  C,  >  1.  all 
the  eigenvalue  trajectories  are  monotonicallv  decreasing  as  R  in¬ 
creases — i.e.,  as  the  inner  load  sphere  increases  in  size.  On  the  other  hand, 
for  ielative  dielectric  strength  er  <  1.  the  eigenvalue  trajectories  increase 
monotonicallv  as  the  inner  sphere  grows.  Increasing  the  index  p  by  unity 
adds  a  knee  to  the  TE  as  well  as  TM  trajectories  for  cr  >  I.  Although  for 
large  er  the  TM  trajectories  do  not  always  show  this  knee,  it  is  present  for 
small  enough  n.  For  er  <  1.  an  identation  is  added  instead. 

Note  that  the  knee  and  the  indentation  shift  the  trajectory  toward  lower  ka 
values  for  TE  and  TM  eigenvalues,  but  preserve  the  monotonicity  of  the 
trajectory. 

We  have  seen  that  regardless  of  the  size  of  the  internal  sphere  loading  the 
dielectric-filled  spherical  cavity,  and  for  the  chosen  values  of  the  dielectric 
strengths  of  the  load  sphere  and  the  medium  filling  the  cavity,  the  lowest- 
lying  mode,  or  fundamental  mode,  of  the  system  will  be  the  TM  eigen¬ 
value  k\]]  (R:£r)a. 

Some  quite  remarkable  behavior  is  manifested  by  the  eigenvalues  as  the 
size  of  the  inner  dielectric  sphere  loading  the  cavity  varies.  We  observed 
short  and  long  ranges  in  the  size  of  the  inner  dielectric  sphere  for  which 
the  eigenvalue  exhibits  virtually  no  significant  change  at  all,  as  well  as 
such  ranges  of  R  over  which  the  eigenvalue  changed  at  a  very  rapid  rate. 
For  many  of  the  higher  order  eigenvalues  several  such  regions  of  behavior 
have  been  seen  following  one  after  another.  This  is  an  unanticipated 
characteristic  for  the  cavity  system. 


It  is  of  interest  to  compare  the  rich  variety  of  behavior  seen  for  the  eigen¬ 
values  of  the  dielectric  loaded  cavity  with  the  homogeneously  filled  con¬ 
centric  conducting  cavity  eigenvalues  1 1 1.  Clearly  for  small  values  of  R. 
as  is  to  be  expected,  the  two  different  cavity  systems  show  the  same  kind 
of  behavior.  Except  for  lower  order  TM  eigenvalues  a  with  er  >  1,  the 
eigenvalues  for  the  two  cavity  systems  are  markedly  different  in  their  de¬ 
pendence  on  R,  the  size  of  the  inner  interface  radius. 

Up  to  this  point,  we  have  constrained  all  our  considerations  to  individual 
families  of  eigenvalue  trajectories  for  the  different  er  values:  i.e.,  we  have 
fixed  the  index  p  and  examined  the  set  of  curves  for  various  ranges  of  the 
index  n.  One  very  clear  feature  of  these  sets  is  that  the  curves  in  each  fam¬ 
ily  of  trajectories  do  not  intersect.  Unfortunately,  this  does  not  yet  give  a 
complete  enough  picture  of  the  eigenvalue  trajectory  characteristics.  We 
shall  now  proceed  to  consider  the  features  of  the  eigenvalues  and  their 
variations  as  R  and  er  change  by  superimposing,  for  each  given  er,  first  all 
the  TE  and  then  all  the  TM  eigenvalues. 

Overview  of  the  TE  and  TM  Eigenvalues  for  er  =  3.00 
We  have  redrawn  in  figure  27  a  superposition  of  the  previously  given  17 
TE  eigenvalue  trajectories  for  the  inner  sphere  of  dielectric  constant  e2  = 
3e,.  It  now  is  immediately  evident  that  many  of  the  trajectories  do  inter¬ 
sect  one  another.  It  is  very  important  to  emphasize  that  each  point  of  inter¬ 
section  of  any  two  eigenvalue  trajectories  represents  a  degenerate  pair  of 
eigenvalues.  And  this  degeneracy  occurs  for  the  given  values  of  Er  at  spe¬ 
cific  sizes,  R ,  of  the  internal  load  sphere.  Figure  27  shows  explicitly  only 
TE-TE  intersections. 

Let  us  begin  a  systematic  examination  of  the  patterns  associated  with  the 
intersection  points  of  the  eigenvalue  trajectories  for  the  k^  (R  x,  =3 )a. 
We  first  note  that  for  er  =  3,  k\\*  a  and  k^a  are  not  degenerate  with  any 
other  TE  eigenvalues.  However,  k2y  a  and  k\[2  a  approach  each  other  at  the 
single  knee  of  k2J  a,  and  the  bottom  of  the  second  slow  decrease  with  ris¬ 
ing  R  in  k\ y<i  with  increasing  R.  In  fact,  we  observe  this  as  a  general  pat¬ 
tern  with  k^a  lying  close  to  k'J\2a  in  the  corresponding  region.  As  n  in¬ 
creases.  this  separation  decreases.  Thus,  at  n  =  5  the  two  trajectories  al¬ 
most  intersect  one  another.  Although  it  is  not  explicitly  shown  for  n  >  5,  a 
pair  of  intersections  will  indeed  occur  in  this  corresponding  region  of  R 
and.  as  n  increases,  the  pair  of  intersections  will  move  apart.  This  same 
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for  the  pairs  k^a  and  o  starting,  for  er  =  3,  with  n  =  2  and  for  all  p. 

A  similar  characteristic  also  occurs  for  the  trajectories  of  k(2V  a  and  k\^  a  at 
the  second  knee  and  the  region  of  R  corresponding  to  the  third  slow 
decrease  as  R  increases.  This  w  ill  be  true  in  general  for  the  pairs  of  trajec¬ 
tories  k'^a  and  a  for  n  >  2.  Clearly  for  the  higher  order  trajectories 

this  pattern  of  behavior  w  ill  be  repeated  at  larger  R  values  corresponding 
to  higher  lying  knees  on  k'^o  and  higher  lying,  slowly  varying  regions  of 

Next  consider  the  eigenvalue  trajectory  pair  of  k^a  and  k\\\i.  It  is  plain 
to  see  in  figure  27  that  they  are  degenerate  at  two  quite  distinct  values  of 
R.  What  we  are  seeing  is  that  there  are  two  sizes  for  the  inner  dielectric 
sphere  at  which  the  system  has  k^a  =k{^a.  One  is  at  R  just  under  0.6 
and  the  other  is  at  R  just  under  0.7.  After  a  little  study  of  figure  27.  we  ob¬ 
serve  that,  for  n  >  1,  the  trajectories  for  k^a  and  k'^2j>-  \ (l  are  similarly 
degenerate  at  higher  R  values,  in  pairs  of  R.  Clearly  these  eigenvalue 
trajectories  are  degenerate  at  two  and  only  two  R  values. 
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Figure  27.  A  composite  of  the  TE  eigenvalues  k^(R: e,  =  })<i  for  the  spherical  cavity  of  radius  a  filled  with 
dielectric  E,  internally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  E2  =  JE,. 
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It  should  be  noted  that  the  trajectories  for  k^a  and  k\\*  a  intersect  at  only 
the  two  points  already  indicated  and  the  trajectories  have  no  other 
degeneracies.  The  trajectories  for  k^a  and  a  behave  in  the  very  same 
manner.  From  this  point  on,  the  trajectories  of  higher  order  radius  will 
display  a  larger  number  of  degeneracies.  Thus  the  trajectory  for  k1^ a  in¬ 
tersects  that  for  *‘V«  twice  and  that  for  k\]]a  twice.  The  k\\  a  trajectory 
also  intersects  the  k^a  trajectory  four  times.  This  degree  of  degeneracy 
becomes  high  rather  rapidly.  Thus  for  example,  the  trajectory  for  k'^o  is 
degenerate  with  the  k^a  trajectory  for  six  different  sizes  of  the  inner 
dielectric  load  sphere.  It  should  be  evident  that  as  ka  exceeds  about  ka  ~  8 
for  mid-range  values  of  R.  the  density  of  TE-TE  degenerate  modes  can  be¬ 
come  rather  large. 

Before  leaving  the  TE  eigenvalues  for  Er  =  3,  we  point  out  a  very  impor¬ 
tant  feature.  Although  the  characteristics  of  the  eigenvalue  trajectories 
may  become  quite  interesting  in  detail  for  0  <  R  <  !,  the  sequence  of 
eigenvalues  at  R  =  1  is  exactly  the  same  as  the  sequence  in  which  the 
eigenvalues  occur  at  R  =  0.  In  general  there  is  a  downward  shift  in  the 
eigenvalues  as  R  increases  from  zero  up  to  unity. 

Next  let  us  examine  the  twenty-one  TM  eigenvalue  trajectories  which 
were  expounded  for  er  =  3  in  about  the  same  detail.  We  first  observe  that 
the  eigenvalue  trajectories  k\\ya  and  k^a  exhibit  no  degeneracy  whatever; 
k^a  has  only  two  degenerate  values,  both  with  kfj  (l ■  Note  that  k1^  a  also 
has  two  degenerate  values,  both  with  k^a.  But  k'fja  also  has  two 
degenerate  values  with  k^ti  for  a  total  of  four  degeneracies  (at  values  of  R 
of  about  0. 1,  0.4,  0.8,  and  0.99).  The  trajectory  for  kfJ a  exhibits  four 
degeneracies  and  that  for  kQa  also  shows  four  sizes  of  internal  sphere  at 
which  degenerate  eigenvalues  occur.  As  we  go  to  higher  order  indices  n,p 
the  number  of  degeneracies  rises  rather  rapidly.  Thus  a  has  two 
degeneracies  with  each  of  the  trajectories  for  k^  a,k^a.  and  k^a  for  a 
total  of  six  degenerate  sets. 

Now  we  observe  that  for  values  of  ka  >  7  the  density  of  TM-TM 
degeneracies  can  become  rather  large.  This  holds  for  nearly  the  entire 
range  of  R.  Again  we  note  that  the  sequence  of  TM  eigenvalues  at  R  =  0  is 
preserved  at  R  =  1  regardless  of  the  behavior  over  intermediate  R  values. 
A  general  shift  to  lower  ka  values  is  evident  for  the  eigenvalues  as  R 


increases. 


We  have  thus  far  only  separately  considered  TE-TE  and  TM-TM 
degeneracies.  The  complete  picture  requires  taking  into  account  TE-TM 
degeneracies  for  er  =  3.  We  now  proceed  to  consider  this  aspect  of  the 
eigenvalue  distribution.  We  will  carry  this  out  without  attempting  to  show 
a  superposition  of  all  the  TE  and  TM  eigenvalues  in  one  figure. 

Now  the  only  eigenvalue  trajectory  that  does  not  have  any  degeneracies  is 
that  for  the  TM  a.  The  TE  eigenvalue  k\\'  a  and  the  TM  eigenvalue 
jfcij1  a  are  both  degenerate,  intersecting  at  about  R  =  0.45  and  at  about  R  = 
0.8.  These  are  the  only  degeneracies  for  these  eigenvalues.  These  two 
eigenvalue  trajectories  are  the  next  lowest  lying  ones,  with  a  lying  be¬ 
low  k\ \] a  for  0  <  R  <  0.45  and  0.8  <  R  <  I  and  interchanging  their  order 
for  the  intermediate  range  of  R.  The  trajectory  for  the  TM  eigenvalue  k^  a 
is  not  only  doubly  degenerate  with  the  TM  eigenvalue  k^a,  but  also  inter¬ 
sects  the  TE  eigenvalue  trajectory  k^a  twice,  first  at  about  R  =  0.5  and 
then  at  about  R  =  0.85.  The  order  of  occurrence  of  these  eigenvalues  at  R 
=  0  is  preserved  at  R  =  1.  In  the  intermediate  range  of  R.  examination  of 
figures  27  and  28  shows  that  their  order  changes  as  R  increases.  Inspec¬ 
tion  of  figures  27  and  28  shows  that  once  we  get  beyond  the  TM  eigen¬ 
values  k1^ a  and  k^u  and  the  TE  eigenvalue  it* '/a  a  considerable  amount 
of  degeneracy  appears,  along  with  shifting  around  in  sequence  of  the 
higher  order  eigenvalues.  Furthermore,  as  ka  increases  into  this  region, 
there  appears  a  rapid  rise  in  the  density  of  the  combined  TE  and  TM  eigen¬ 
values.  Careful  inspection  of  the  figures  reveals  a  considerable  regularity 
in  the  occurrence  of  intersections  of  the  set  of  combined  trajectories.  We 
will  not  elaborate  any  further  on  this  subject. 

Having  seen  the  considerable  complexity  in  the  set  of  combined  TE  and 
TM  eigenvalue  trajectories  for  er  =  3.  we  next  turn  to  considering  the 
changes  that  result  from  varying  the  strength  of  the  dielectric  sphere  load¬ 
ing  the  cavity. 
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THE  m  EIGENVALUES  k£'  (R 

Figure  28.  A  composite  of  the  TIM  eigenvalues  kl^(R\?r  =3)a  for  the  spherical  cavity  of  radius  a  filled  with 
dielectric  E,  internally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  £2  =  JE,. 

5.3  Overview  of  the  TE  and  TM  Eigenvalues  for  er  =  10.00 

Next  we  shall  examine  the  superposition  in  figure  29  of  the  set  of  twenty- 
two  TE  eigenvalues,  which  we  considered  individually  earlier.  We  now 
assume  the  inner  dielectric  sphere  to  be  of  strength  e2  =  l()e,.  In  figure  29 
for  er  =  10,  we  see  that  only  the  eigenvalue  trajectory  k\\ a  has  no  intersec¬ 
tions  with  any  other  TE  trajectory.  Recall  that  for  er  =  3  there  were  two 
such  eigenvalue  trajectories. 

The  kW'iRXr  =  10)  cf  trajectory  nearly  coincides  with  the  k\\  (R:er  =  3)  a 
trajectory  for  the  approximate  range  0  <  R  <  0. 15.  For  higher  R  values,  the 
k\\ (R:£r  =  10) a  trajectory  lies  lower.  We  also  note  that  the  eigenvalue 
trajectory  for  k\\’  a  at  its  knee  lies  near  the  eigenvalue  trajectory  for  k^a 
just  beyond  its  first  knee  and  where  it  begins  its  slowly  decreasing  region 
with  increasing  R.  This  is  true  in  the  same  R  region  in  general  for  the 
eigenvalue  pairs  k)}j  a  and  for  all  n  >  1.  Although  we  have  not  ex¬ 
plicitly  shown  it  in  figure  29,  there  will  be  a  value  of  n  where  this  set  of 
trajectories  intersects  and  continues  to  do  so  for  all  higher  values  of  the  in¬ 
dex  n.  Many  similar  pairs  of  trajectories  exhibit  this  same  behavior,  as  can 
be  seen  by  careful  study  of  figure  29.  We  shall  not  go  into  this  in  any  fur- 
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ther  detail.  We  shall  instead  examine  a  further  complexity  in  the 
degeneracy  of  the  eigenvalues. 


The  TE  eigenvalue  trajectory  k^a  intersects  the  TE  trajectory  for  k'^a  for 
two  sizes  of  internal  dielectric  sphere  load — at  about  R  =  0.15  and  at  about 
R  -  0.25.  Note  that  k^u  shows  TE-TE  degeneracy  only  at  these  two 
values.  The  eigenvalue  trajectory  for  k\\' a  has  two  points  of  degeneracy 
with  the  trajectory  for  k\^ a — at  about  R  =  0. 13  and  R  =  0.35.  It  also  has 
two  points  of  degeneracy  with  the  trajectory  for  a — at  about  R  =  0.2 
and  R  =  0.3.  Thus  the  eigenvalues  k\\]  a  and  k^a  each  have  TE-TE 
degeneracy  at  four  sizes  of  internal  sphere  load.  Note  that  the  points  of 
degeneracy  for  these  occur  only  for  relatively  small  internal  spheres.  The 
situation  really  increases  in  complexity  very  rapidly  as  we  go  to  higher  or¬ 
der  eigenvalues.  Thus  the  eigenvalue  k^'a  displays  six  TE-TE  degenerate 
situations,  and  k^u  also  shows  six  such  degeneracies.  The  eigenvalue 
k<\  a  shows  10  TE-TE  degeneracies.  Clearly  the  number  of  TE-TE  degen¬ 
eracies  is  growing  rapidly.  We  note  that  up  to  about  ka  =  7,  the  density  of 
TE-TE  degeneracies  is  rather  substantial,  but  only  for  R  values  up  to  about 
0.5.  For  higher  ka  values,  this  high  density  of  TE-TE  degeneracies 
spreads  to  much  larger  R  values.  What  we  are  seeing  is  the  general 
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tendency  of  the  eigenvalues  to  move  to  lower  ka  as  er  increases,  with  the 
higher  order  trajectories  moving  farther. 

Figure  30  contains  the  superposition  of  all  the  twenty-three  TM  eigenvalue 
trajectories,  k^(Rxr  =  10)  tf.  which  we  examined  earlier  in  sets  of  eigen¬ 
values  for  fixed  index  p.  Let  us  examine  figure  30  for  TM-TM 
degeneracy.  The  First  feature  to  note  is  that  only  one  eigenvalue  trajectory 
does  not  intersect  any  other,  namely,  the  k\\'(R:£r  =  10) a.  Thus,  in  going 
from  er  =  3  to  er  =  10,  the  number  of  eigenvalues  that  have  no  trajectories 
with  TM-TM  degeneracy  is  reduced  from  two  to  only  one.  The  eigen¬ 
value  curve  k{2\  a,  which  exhibits  no  TM-TM  degeneracy  for  er  =  3.  now 
has,  for  er  =  10,  degeneracy  with  A(|V  ( Rxr  =  10)  at  about  R  =  0.35  and  R 
=  0.65.  These  intersections  are  the  only  internal  sphere  radii  at  which  k\~2  a 
has  TM-TM  degeneracy.  There  are  four  TM-TM  degeneracies  for  A,j  a. 
Two  are  with  k^a  near  R  =  0.25  and  R  =  0.95,  corresponding,  respec¬ 
tively.  to  a  small  and  a  large  load  sphere  in  the  interior,  and  two  with  k^u 
are  at  about  R  =  0.30  and  R  =  0.55,  corresponding  in  both  instances  to 
intermediate-size  inner  dielectric  spheres.  Next  we  see  that  the  trajectory 
of  eigenvalues  k$  a  has  six  points  of  intersection,  with  a  total  of  three 
other  TM  eigenvalue  trajectories.  (There  are  two  degeneracies  with  each 
of  k^a,  A^’rt,  and  k™a.)  Also,  the  degeneracies  occur  for  inner  load 
spheres  whose  sizes  spread  over  almost  the  whole  range  of  R.  It  is  of  in¬ 
terest  also  that  for  small  inner  spheres  up  to  about  R  =  0.10.  the  eigen¬ 
values  A*,2,1  a  and  A*,2’  are  very  close  to  one  another. 

Going  up  in  order,  the  eigenvalue  k'^  a  has  a  trajectory  with  eight  intersec¬ 
tions  with  other  TM  eigenvalue  trajectories.  It  has  two  intersections  with 
each  of  k\2l  a,  k\2Ja,  k^a,  and  k^a.  Again,  we  can  readily  see  that  these 
degeneracies  correspond  to  inner  sphere  radii  almost  over  the  whole  range 
of  R.  From  about  R  =  0.60  to  about  R  =  0.70,  k™  a  and  A1,2'  a  are  very 
nearly  equal.  Thus  we  see  the  feature  of  two  ranges  in  R  for  which  the 
eigenvalue  trajectory  for  A*,2,’ a  is  close  to  degeneracy.  If  we  continue  to 
examine  the  higher  order  eigenvalue  trajectories,  e  g.,  A'f/ a.  we  find 
twelve  TM-TM  degeneracies  of  a  large  variety  of  characteristics.  Clearly 
the  TM-TM  degeneracy  picture  is  becoming  considerably  involved  rather 
rapidly  now. 
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Figure  30.  A  composite  of  the  I'M  eigenvalues  k(^(R:tr  -  10) <i  for  the  spherical  cavity  of  radius  a  filled  with 
dielectric  E,  internally  loaded  with  a  concentric  dielectric  sphere  of  radius  h  and  strength  f  2  =  lllf ,. 


Several  novel  features  are  apparent  in  figure  30.  Degeneracy  is  i.uln-i 
common  for  inner  dielectric  spheres  that  nearly  fill  the  entire  (..nits  I: 
also  occurs  frequently  for  quite  small  inner  spheres  that  load  the  ■. 

These  features  were  already  present  for  Er  =  3.  but  are  more  emptv. 
for  the  higher  er  case.  There  are  substantial  ranges  of  R  tor  who.  t  m  , 
TM-TM  degeneracy  occurs.  The  density  of  intersections  of  IM  ttaiev 
tories  becomes  quite  high  rather  soon.  Many  other  regular  features  <.  an  be 
discerned  in  figure  30  which  we  shall  not  explicitly  elaborate  on  Finally 
we  again  note  that  the  order  of  the  eigenvalues  near  R  =  0  is  preserv  ed  in 
the  limit  as  R  goes  to  unity.  The  general  property  of  a  downward  shift  of 
all  trajectories  to  lower  ka  with  increasing  Er  is  evident  at  R  =  0. 
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To  get  a  much  more  thorough  picture  of  the  complicated  situation  that  oc¬ 
curs  at  £r  =  10  really  requires  that  a  single  figure  containing  all  forty-five 
TE  and  TM  eigenvalue  trajectories  superimposed  be  available  for  a  de¬ 
tailed  analysis.  Obviously  such  a  picture  would  be  a  highly  confused  one. 
losing  clarity  and  accuracy  of  detail  rather  quickly.  What  we  shall  use  as 
an  alternative  mechanism  to  acquire  the  more  complete  overview  of  the 
whole  TE  and  TM  situation  is  simultaneous  reference  to  both  figure  29  for 
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the  TE  eigenvalues  and  figure  30  for  the  TM  eigenvalues,  for  £r  =  10.  In 
this  manner,  we  can  study  all  TE-TE,  T'E-TM,  and  TM-TM  degeneracies 
at  once. 

At  the  outset,  we  note  one  very  striking  characteristic.  The  TM  eigen¬ 
values  of  ft)2,’  (R  ,er  =  10)  a,  which  are  the  lowest  for  all  /?,  are  the  only 
ones  not  degenerate  at  any  value  of  R  with  any  other  eigenvalue.  This 
property  is  not  new  by  any  means.  What  does  arise  now  is  the  feature  that 
over  the  range  of  R  from  about  R  =  0.43  to  about  R  =  0.65  the  trajectory 
for  ft)2’ a  lies  very  close  to  the  TE  eigenvalue  trajectory  for  k\\]  a.  This 
raises  an  interesting  and  important  question.  Namely,  if  we  take  e2  suffi¬ 
ciently  larger  than  £,,  can  we  eventually  attain  degeneracy  even  for  what 
so  far  has  been  the  lowest  eigenvalue,  i.e.,  ft)2’ al  Indeed,  it  can  be  demon¬ 
strated  that  for  er  =  100,  for  example,  the  k\\}  a  curve  intersects  the  ft)2,1  a 
curve  at  two  values  of  R  (at  approximately  R  =  0.15  and  R  =  0.92).  Thus, 
for  this  case,  all  the  modes  are  degenerate  for  some  size  of  inner  dielectric 
sphere  that  we  load  the  cavity  with.  Furthermore,  when  £r  =  100,  the 
TM,,  mode  is  the  dominant  mode  only,  for  small  and  large  values  of  R  = 
b/a.  In  the  intermediate  range  of  R ,  which  is  bounded  by  the  two 
TE,,-TM,,  degeneracies,  the  TE,,  mode  surplants  the  TM,,  as  the  domi¬ 
nant  mode! 

It  can  be  further  inferred  that  there  exists  a  value  of  £f  between  10  and  100, 
at  which  the  ft)1/ a  trajectory  just  grazes  the  ft1,2/ a  trajectory,  thus  producing 
only  one  value  of  degeneracy,  and  becoming  the  bounding  Er  below  which 
only  the  TM,,  mode  can  be  the  dominant  mode,  independent  of  R.  If  one 
wishes  to  generalize  further,  "grazing"  values  of  Er  should  exist  for  the 
other  degeneracy  intersections  as  well.  It  would  not  be  a  difficult  exercise 
to  solve  for  these  "grazing"  values  of  £r,  although  we  have  not  addressed 
that  particular  study  for  this  report. 


Continuing  to  investigate  the  situation  for  er  =  10  with  regard  to 
degeneracy  as  we  have  earlier,  we  find  that  the  TE  eigenvalue  curve  ft,1,1  a 
is  degenerate  only  with  the  TM  curve  ft)2’  a  near  R  =  0.2  and  near  R  =  0.95, 
i.e.,  for  a  small  and  for  a  quite  large  inner  sphere.  The  TM  eigenvalue 
curve  k{2\  a  is  degenerate  with  other  modes  for  six  distinct  inner  sphere 
sizes.  In  addition  to  the  degeneracies  with  the  TE  mode  ft/,’ a,  it  is  twice 
degenerate  with  the  TE  mode  ft/,’ a  and  also  with  the  TM  mode  ft)2’  a.  Note 
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that  the  lower  R  values  (about  R  =  0.35)  at  which  the  a  and  occur 
are  nearly  equal. 

Proceeding  to  higher  order,  the  TM  eigenvalue  curve  k^' a  has 
degeneracies  at  10  values  of  R.  In  addition  the  trajectory  for  k\ j'</  nearly 
coincides  w  ith  several  other  eigenvalue  trajectories  for  several  ranges  of  R. 
Thus  once  we  get  beyond  about  R  =  0.1,  the  density  of  degeneracies  grows 
very  rapidly.  Further,  as  we  go  to  higher  order  eigenvalues,  this  lower 
threshold  in  R  decreases  steadily.  Clearly  for  er  =  10  we  find  a  very  dense 
distribution  of  degeneracies  and  a  rich  variety  of  near  degeneracy  as  well. 
Again,  the  sequence  of  TE  and  TM  eigenvalues  near  R  -  0  remains  in¬ 
variant  as  R  approaches  the  limit  of  unity,  although  many  interesting  varia¬ 
tions  in  behavior  occur  in  between  these  limits.  Recall  also  that  all  the  TE 
and  TM  eigenvalues  decrease  monotonically  as  R  increases  from  zero  to 
unity,  where  the  inner  sphere  fills  the  entire  cavity.  The  general 
downward  shift  to  lower  ka,  for  all  the  trajectories,  which  is  larger  for 
larger  er,  tends  to  increase  the  density  of  degeneracies  for  R  >  0. 

We  have  now  accumulated  a  considerable  body  of  knowledge  concerning 
the  degeneracies  in  our  system  and  the  effect  of  increasing  er  on  them. 
Next  we  consider  the  corresponding  situation,  where  we  fill  the  cavity 
with  a  high  dielectric  medium  and  then  load  it  in  its  interior  with  a  con¬ 
centric  dielectric  sphere  of  lower  dielectric  strength,  i.e.,  er  <  1.  We  shall 
examine  the  results  only  in  the  single  case  of  er  =  1/3. 

Overview  of  the  TE  and  TM  Eigenvalues  for  er  =  1/3 
We  shall  begin  as  before,  by  first  examining  the  TE-TE  degeneracy  situa¬ 
tion.  Figure  3 1  contains  in  superposition  all  twenty  of  the  TE  eigenvalue 
trajectories  we  examined  as  single  sets  of  curves  for  fixed  values  of  the  in¬ 
dex  p.  Unlike  what  we  observed  in  figures  27  and  29,  we  now  see  that  the 
first  three  TE  trajectories,  k\\ a,  k^  a,  and  k\\  a,  are  free  of  degeneracy  for 
all  R.  Furthermore,  these  eigenvalues  are  well  separated  from  all  other  TE 
eigenvalues  for  the  entire  range  of  R.  The  k'^a  curve  shows  the  first 
degeneracy.  Its  trajectory  intersects  the  k\\' a  trajectory  at  about  R  =  0.5 
and  again  at  about  R  =  0.85.  The  eigenvalues  of  k\\ a  are  degenerate  at 
four  values  of  R — twice  with  k\^a  and  twice  with  k^a.  All  the  TE-TE 
degeneracies  cited  thus  far  are  well  separated  in  R  from  each  other.  The 
trajectory  for  the  eigenvalues  k\^  a  shows  six  sizes  of  inner  sphere  at 
which  degeneracy  occurs.  There  are  two  values  of  R  for  each  of 
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leu  a,  a,  and  a  at  which  these  intersections  lie.  The  k^u  trajectory 
also  has  six  degenerate  values  of  R ,  two  with  each  of  the  trajectories  for 
$  a  and  k^  a,  and  although  it  is  not  explicitly  shown,  two  more  with 
k\\  a.  The  a  trajectory  intersects  k^ a  twice  and  k\ '2'  a  twice. 

The  k^\  a  and  kj^ci  trajectories  are  nearly  coincident  for  /?,  ranging  all  the 
way  up  to  beyond  R  =  0.3.  Similar  behavior  exists  for  k\\]  a  and  k^o. 
There  are  six  degenerate  R  values  for  k^\  a,  two  with  each  of  k\'j  a ,  kjj  cu 
and  k\j  ci.  The  next  trajectory,  k\'J  a,  has  eight  values  of  R  at  which 
degeneracy  occurs.  We  can  observe  that  the  density  of  degenerate  pairs  of 
eigenvalues  increases  with  increasing  orders  of  n  and  p ,  but  at  a  slower 
rate  than  we  saw  happen  for  £r  >  1.  Indeed,  up  to  ka  ~  6.5,  very  little 
degeneracy  can  be  found.  Inspection  of  figure  31  reveals  that  as  R  in¬ 
creases  the  eigenvalues  are  shifted,  in  general,  to  higher  values.  Again  we 
note  that  the  sequence  of  eigenvalues  at  R  =  0  is  preserved  at  R  =1. 


Next  we  examine  the  thirty-one  TM  eigenvalue  trajectories  for  £r  =  1/3  all 
superimposed  in  figure  32.  We  begin  by  noting  that  the  four  lowest  order 
TM  eigenvalues  of  kff  a,  k^a,  kfj a,  and  k^ a  are  not  degenerate  for  any 


Figure  31.  A  composite  of  the  TK  eigenvalues  k^(R,er  =  1/3) «  for  the  spherical  cavity  of  radius  a  filled  with 
dielectric  £,  internally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  E2  =  !/3e(. 


values  of  /?.  Further,  there  is  a  clear  separation  of  all  these  eigenvalues  for 
all  R  values,  with  several  exceptions.  One  exception  is  for  R  near  unity, 
where  k™ a  is  very  close  in  value  to  kf]a.  The  second  exception  is  for  the 
range  of  R  up  to  about  R  =  0.15,  where  k^a  and  k{2]  a  nearly  coincide  in 
value.  The  TM  eigenvalues  k^a  and  kfl  ci  are  degenerate  with  each  other 
at  two  R  values,  the  lower  at  about  R  =  0.4  and  the  upper  just  under  R  = 
0.9.  The  eigenvalue  trajectory  for  k^  a  also  shows  two  inner  sphere  sizes 
where  degeneracy  occurs.  These  are  at  the  values  of  R  about  R  =  0.45  and 
at  about  R  =  0.95.  Both  of  these  intersections  are  with  the  trajectory  for 
the  eigenvalue  k^a.  In  addition  to  these  two  degenerate  R  values  with 
k^’  a.  the  eigenvalues  k^a  have  two  degenerate  values  of  R  with  k^a. 
These  latter  R  values  are  at  about  R  =  0.55  and  R  =  0.95.  The  values  at 
which  the  degeneracies  occur  for  each  eigenvalue  are  well  separated  in  R. 

The  trajectory  for  k^a  for  six  values  of  R  is  degenerate  with  other  eigen¬ 
values.  These  include  two  trajectory  intersections  with  k\2^  a.  One  is  at 
about  R  =  0.05  and  the  other  very  close  to  R  =  1.0.  Two  more  are  with 
kf-^cJ — one  near/?  =  0.4  and  the  second  near/?  =  0.95.  The  remaining  two 
values  are  the  degeneracies  with  kj^a  already  indicated.  The  eigenvalue 


Figure  32.  A  composite  of  the  TM  eigenvalues  k(^(R:er  =  1/3)  a  for  the  spherical  cavity  of  radius  a  Tilled  with 
dielectric  E,  internally  loaded  with  a  concentric  dielectric  sphere  of  radius  b  and  strength  E2  =  1/36,. 
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trajectory  k^ti  also  has  six  sizes  of  inner  sphere  at  which  degeneracy  ex- 
sists  with  other  eigenvalues.  Two  with  a  have  already  been  pointed 
out.  Two  more  occur  with  each  of  two  eigenvalues  from  k^]’  a  and 
These  are  well  spread  out  over  R  values  from  0.5  up  to  0.9.  All  higher  or¬ 
der  eigenvalue  trajectories  will  show  increasing  numbers  of  R  values  at 
which  degeneracy  occurs  with  a  larger  number  of  other  eigenvalues.  Ex¬ 
amination  of  figure  32  shows  that  above  ka  =  6  the  density  of  degeneracies 
grows  rather  rapidly.  Furthermore,  we  see  in  figure  32  the  increase  in  fre¬ 
quency  of  occurrence  of  ranges  of  R  over  which  many  pairs  of  eigenvalues 
nearly  coincide.  These  ranges  are  quite  large  in  many  cases. 

As  we  expec‘.  the  order  of  the  TM  eigenvalues  at  R  =  0  is  preserved  at  R  = 
1 .  There  is  present  a  genera!  upward  shift  in  the  eigenvalue  trajectories  for 
R  >  0  for  £r  <  I .  We  recall  that  now  all  the  eigenvalues  are  monotonically 
increasing  as  R  increases.  Many  other  interesting  features  are  present  in 
figure  32  showing  all  the  TM  trajectories  for  er  =  1/3.  We  shall  not  go  any 
further  into  this  subject,  but  shall  instead  go  on  to  consider  the  composite 
of  all  fifty-one  eigenvalue  trajectories  for  Er  =  1/3. 

For  obvious  reasons,  we  do  not  include  such  a  superposition  plot  of  all  the 
eigenvalues.  We  can  analyze  all  the  eigenvalue  trajectories  by  simul¬ 
taneously  examining  figures  31  and  32. 

The  very  first  observation  to  be  made  is  that  the  TM  trajectories  k^  a  and 
k1^  a  are  nowhere  degenerate  with  any  other  eigenvalues  for  any  R  value. 
For  all  R,  the  eigenvalues  of  k{^a  and  k£\ a  are  well  separated  from  all 
other  eigenvalues.  The  TE  eigenvalues  of  k\^  a  and  TM  eigenvalues  of 
£3,  a,  which  respectively  had  no  TE-TE  and  no  TM-TM  degeneracy,  inter¬ 
sect.  There  is  TE-TM  degeneracy  of  these  two  curves  at  about  R  =  0.55 
and  at  about  R  =  0.9.  Now  we  note  that  the  TE  curve  k^a  has  two  values 
of  R  at  which  it  is  degenerate  with  the  TM  curve  k™  a.  These  values  near 
R  =  0.5  and  near  R  =  0.97  are  the  only  degeneracies  in  k^  a.  The  TM 
eigenvalue  k4 ,  a  is  also  degenerate  at  only  two  values  of  /?,  near  R  =  0.5 
and  R  =  0.97.  Both  degeneracies  are  with  the  TE  eigenvalue  curve  k^a. 
As  we  have  already  noted,  k^  a  is  nearly  coincident  with  k^  a  for  R  up  to 
about  R  =  0.15  and  again  near  R  =  1 . 

There  are  four  R  values  at  which  the  TM  eigenvalues  of  kfj  a  are 
degenerate.  Two  are  with  the  TE  eigenvalues  of  k\\]  a  at  about  R  =  0.35 


and  also  at  about  R  =  0.7.  and  two  are  with  the  TM  eigenvalues  of  a  at 
about  R  =  0.4  and  about  R  =  0.9.  The  k{3\]  a  trajectory  is  degenerate  at  four 
R  values,  the  two  already  referred  to  with  k^u  and  two  others  with  the 
eigenvalues  of  k\]'  a  at  about  R  =  0.5  and  R  =  0.98.  Four  degeneracies  ex¬ 
ist  for  k[]'  a,  the  pair  already  mentioned  with  the  TM  curve  k'flu  and  the 
two  with  the  TE  curve  k\\' a.  From  this  point  on,  the  higher  order  eigen¬ 
values  show  rapid  increase  in  degeneracy.  The  density  of  degeneracy  very 
quickly  becomes  high  beyond  ka  =  4.  In  addition,  many  eigenvalue  trajec¬ 
tories  are  nearly  coincident  with  other  eigenvalue  trajectories  for  various 
ranges  of  R.  Some  ranges  of  R  are  quite  extensive  where  this  situation  oc¬ 
curs.  Thus,  for  example,  the  TE  curve  k\\'  a  and  the  TM  curve  k\^  a  are 
nearly  coincident  for  values  of  R  from  zero  up  to  about  R  =  0.5. 

All  the  er  =  1/3  trajectories  increase  monotonically  with  R.  Finally,  we 
state  the  expected  condition  that  the  sequence  in  which  the  eigenvalues  are 
ordered  at  R  =  0  remains  invariant  at  R  =  1 .  As  usual,  there  are  many  in¬ 
teresting  features  present  in  the  complete  set  of  eigenvalue  trajectories  for 
er  =  1/3  which  we  shall  not  go  into  explicitly. 

Final  Summary  and  Discussion 

We  have  analyzed  the  results  of  solving  equations  (32)  and  (36)  for  the  TE 
and  TM  eigenvalue  trajectories,  for  a  dielectric-filled  spherical  cavity 
loaded  with  an  inner  concentric  dielectric  sphere. 

By  choosing  the  inner  dielectric  sphere  of  higher  strength  than  that  of  the 
cavity  dielectric  for  various  values  of  Er  >  1  and  reversing  the  roles  of  the 
dielectrics  and  thus  obtaining  er  <  1,  we  studied  the  effects  of  eigenvalue 
trajectories.  The  specific  values  chosen  for  er  were  3.00,  10.1X),  and  1/3. 
These  constitute  a  reasonably  large  range  for  er.  For  values  of  er  >  1  all 
the  eigenvalues  decrease  monotonically  as  the  radius  of  the  inner  sphere 
increases.  On  the  other  hand,  for  er  <  1,  ail  the  eigenvalues  increase 
monotonically  as  the  radius  of  the  inner  sphere  increases.  In  all  cases  of 
er,  although  the  eigenvalues  display  some  very  interesting  behavior  for  in¬ 
ner  sphere  radii  between  R  =  0  and  R  =  1 ,  the  eigenvalue  sequence  of  R  = 
0  is  preserved  at  R  =  1.  Comparing  two  cases  of  er  >  1,  we  find  a  general 
shift  of  the  eigenvalues  to  lower  values,  for  R  >  0,  for  the  higher  er  case 
relative  to  the  lower  er  case.  For  two  values  of  er  <  1 ,  the  reverse  situation 
holds  true.  Varying  er  induces  a  multitude  of  changes  in  which  modes  of 
the  system  are  degenerate  and  for  which  values  of  radius  of  the  inner 
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sphere  the  degeneracy  occurs.  Further  varying  er  affects  the  density  of 
degeneracies  among  the  eigenvalues. 

A  number  of  practical  applications  of  the  results  of  this  study  suggest 
themselves.  For  example,  we  see  explicitly  for  the  simple  system  studied 
that  the  characteristics  of  ait  empty  cavity  and  a  loaded  cavity  can  be  enor¬ 
mously  different.  We  shall  show  in  subsequent  reports  that  the  interior 
contents  of  a  cavity  can  be  sensed  in  the  far  Field.  The  results  reported  on 
in  this  paper  indicate  that,  if  we  appropriately  add  elements  to  the  interior 
of  the  cavity,  identification  of  the  cavity  and  its  interior  components  can  be 
made  very  difficult.  Thus  we  noted  that  adding  a  dielectric  sphere  of  high 
strength  inside  shifts  many  high  order  modes  to  lower  values.  This  gives 
rise  to  confusion  in  the  array  of  resonances  in  the  cavity  rather  than  the 
simple  response  of  only  a  few  resonances.  We  also  see  analytic  proof  of 
the  large  number  of  resonances  that  enter  into  the  cavity  response  when 
the  cavity  contains  even  a  very  simple  element.  The  results  of  this  study 
show  that  it  can  be  very  difficult  to  achieve  and  maintain  EMC  and  EMI 
hardening  of  systems  by  introducing  protective  devices  for  selected  fre¬ 
quencies  because  of  the  large  number  of  frequencies  that  must  be  account¬ 
ed  for  when  even  simple  elements  are  in  the  interior  of  metallic  shields. 
We  have  further  shown  that  the  well  separated  lower  modes  of  systems 
can  be  shifted  to  higher  frequencies  if  we  coat  the  inner  surface  of  the  ex¬ 
terior  metal  boundary  with  high  dielectric  strength  materials.  We  have  by 
no  means  exhausted  all  possible  application  potentials  of  this  study. 

One  final  important  area  of  application  that  suggests  itself  in  this  study  is 
that  of  the  difficult  mode  conversion  problem.  Using  the  work  of  Rowe 
and  Waiters  |3|  and  that  of  Kindermann  |4|,  Doane  |5|  showed  how  the 
degeneracy  of  the  TE01  and  the  TMn  circular  waveguide  modes  could  be 
used  to  convert  from  the  first  to  the  second.  The  high  degree  of  flexibility 
we  have  seen  in  controlling  degeneracy  among  the  many  modes  in  the 
study  we  reported  on  presents  us  with  an  opportunity  to  develop  an  addi¬ 
tional  mode  conversion  capability. 

We  conclude  the  discussion  by  pointing  out  that  one  can  achieve  a  deeper 
understanding  of  what  we  found  in  the  study  reported  on  here  by  studying 
the  behavior  of  the  spatial  distribution  of  the  field  components  in  our  sys¬ 
tem.  This  will  be  presented  in  a  series  of  subsequent  reports. 
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ATTN  RECORD  COPY,  A3NC-LAB-T3 
ATTN  TECHNICAL  REPORTS  BRANCH, 
ASNC-LAB-TR  (3  COPIES) 

HARRY  DIAMOND  LABORATORIES 

ATTN  D/DIVISION  DIRECTORS 

ATTN  LIBRARY,  SLCHD-Tl  (3  COPIES) 

ATTN  HDL  LIBRARY,  SLCHD-TL  (WOODBRIDGE) 

ATTN  ZABLUDOWSKI,  B.,  SLCHD-IT-EB 

ATTN  JOHNSON,  T. ,  SLCHD-ST 

ATTN  PISANE,  G.,  SLCHD-ST-AB  (20  COPIES) 

ATTN  RUTH,  B.,  SLCHD-NW-RI 

ATTN  LIBELO,  L.,  SLCHD-NW-RI  (30  COPIES) 
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